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Abstract: We compute the differential B — > decay width in the Standard Model as a 
function of the photon energy using Dressed Gluon Exponentiation (DGE). The resummed 
spectrum is matched with the fixed-order expansion, making use of the next -to-next-to- 
leading order (NNLO) results for the matrix element of the magnetic dipole interaction O7 
and NLO ones for other operators in the effective Weak Hamiltonian. We develop a new 
technique to implement constraints on the analytic structure of the Sudakov factor in mo- 
ment space. This improves the behavior of the resummed spectrum away from the Sudakov 
region. We also derive an analytic expression for the Borel transform of the perturbative 
series for the 0^ spectrum in the large-/3o limit. Using this example we demonstrate that 
exponentiation in moment space is necessary for the calculation of the spectrum for ^ 2 
GeV. Finally, we investigate numerically the relation between renormalons, power correc- 
tions and support properties. We present predictions for the branching fraction and the 
first few spectral moments as a function of a cut > £'0 and estimate the theoretical 
uncertainty. 
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1. Introduction 

Inclusive radiative B decays, B — > Xgj, have become an essential ingredient in precision 
tests of the Standard Model. The Standard Model 6 — > 57 decay occurs only through 
loops (penguin diagrams) involving the W Boson, whose mass is significantly larger than 
the available energy, mw S> mi,. This makes the B — > Xg^ width a sensitive probe of 
any potential flavor-changing short-distance interaction beyond the Standard Model, see 
e.g. [1-3]. 

The Standard Model Branching Fraction (BF) is known [4-8] since a few years to 
next-to-leading order (NLO) in renormalization-group improved perturbation theory. The 
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summation of large logarithms of mff/mw is conveniently formulated in the framework of 
the effective Weak Hamiltonian, where virtualities of order of the Weak scale are integrated 
out to obtain a set of local operators Oi of dimension 6. The result for the inclusive decay 
width with a minimal photon-energy cut, Ej > Eo, takes the form: 

T{B Xsj,E> Eo) = \VM' (mf (m,))' mf (1.1) 

X ^ Cfifi)Cf{^i)G,,{Eo,fi), 

where rrib is the quark pole mass, Gij {Eq, /v,) are matrix elements of operators in the effective 
Weak Hamiltonian, Gij{Eo,ii) = Xlx, (^| 0|(//) 1X^7) (Xs7| Oj(//) |^) where E^ > Eq, 
and C^^{iJ,) are the corresponding Wilson coefficients. The matrix elements can be com- 
puted in perturbation theory, replacing the B meson by an on-shell b quark and the 
hadronic system Xg by a partonic one, owing to the inclusive sum over the final states. 
This replacement can be justified using the Operator Product Expansion (OPE) so long as 
the photon-energy cut is insignificant. Upon removing the cut, power corrections to the 
partonic calculation can be formally shown to be of 0(A^/m^), and they are numerically 
small, approximately +2.5% of the total BE [7]. The Standard Model BE, which has been 
determined with about ±10% uncertainty [7,8], is found to be in good agreement with 
experimental measurements by CLEO, Belle and BaBar [9-14]. 

The current world average of all experimental data, prepared by the Heavy Flavor 
Averaging Group [15], is B{B — > Xsj,E., > 1.6 GeV) = (355±24±10±3)xl0-6, where the 
errors are: combined statistical and systematic uncertainty, "shape-function" uncertainty 
owing to the extrapolation from the region of measurement E^ > Eq, where Eq > 1.8 GeV, 
to the reference range of Ej > 1.6 GeV, and uncertainty owing to the b ^ dj fraction. The 
"shape-function" uncertainty varies significantly between different theoretical approaches. 
For example, Ref. [16] assigns an error of about ±8% to the extrapolation below = 1-8 
GeV, i.e. three times the size of the error quoted here. 

One of the essential ingredients for improving the precision of this comparison in the 
future is the theoretical calculation of the photon energy spectrum. Owing to irreducible 
background, experimental measurements are limited to the range E^ > 1.8 GeV and they 
can be significantly improved if the requirement on the range of measurement is relaxed, 
for example, to E^ > 2.0 GeV. This, however, requires larger extrapolation that relies on 
the theoretical description of the spectrum. Fortunately, this extrapolation presents very 
little sensitivity to short-distance physics, in sharp contrast with the total width. On the 
other hand, it requires detailed understanding of the QCD dynamics. 

The QCD calculation of inclusive decay spectra is essential for other aspects of flavor 
physics. An important example is the determination of \V-ah\ from inclusive charmless 
semileptonic decays, B — > Xylu, where the background due to the 50 times more abundant 
semileptonic decay into charm restricts the region of measurement to Mx < 1-7 GeV. The 
dynamics there is similar to the one governing the B — > Xg^ spectrum to the extent 
that spectral measurements of the photon-energy spectrum in B — > Xg'y are used for the 
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determination of \V^\^\. In recent years the B — > Xgj spectrum has become the prime 
testing-grounds for theoretical approaches to inclusive distributions, see e.g. [16-29]. 

The main challenge in computing inclusive decay spectra in QCD is the complex dy- 
namics of the threshold region. In this region the decaying b quark is just slightly off its 
mass shell, owing to its "primordial" Fermi motion and to soft gluon radiation [30-33]. 
Specifically, the B — > Xg'j spectrum peaks near the partonic threshold, — > mb/2, or 
X = 2E-y/rrH) — > 1. This region is characterized by parametrically-large higher-order per- 
turbative corrections (Sudakov logs) [33] as well as non-perturbative effects, predominantly 
ones related to the Fermi motion of the b quark in the B meson. 

It is universally acknowledged that fixed-order perturbative results cannot be directly 
used for comparison with spectral data, not even the first few moments of the photon 
energy with experimentally-relevant cuts, such as > 1.8 GeV. Fixed-order results for 
the B — > Xg^y spectrum are characterized by 

• Sudakov logarithms, namely singular real-emission corrections to the differential spec- 
trum dT{x)/dx, of the form hi^(l — x) /(I ~ x) with < 2n — 1 at order a^", owing to 
multiple soft and collinear radiation. The perturbative spectriim is nevertheless intc- 
grable as there are also infrared-singular virtual corrections proportional to 5{1 — x) 
— the spectral moments, 

^"io r-, dx ^'-'^ 

are infrared safe. 

• Support for E^ < mh/2, where is the quark pole mass, setting the upper limit of 
integration in Eq. (1.2) as a; = 1. The perturbative support at any order is different 
from the physical one, E^ < Mb/2, where Mb is the meson mass. Importantly, the 
pole mass itself has a linear infrared renormalon ambiguity [34-36], mi,±0{A), 
and therefore it cannot be assigned a precise value without specifying an additional 
regularization prescription. In a fixed-order framework one computes the pole mass 
order- by-order from a given short-distance mass, such as m^^, but the result strongly 
depends on the order, as the series is badly divergent. The use of alternative mass 
schemes [26, 27, 37, 38] amounts to introducing an infrared cutoff, which hinders the 
possibility of using of the inherent infrared safety of the on-shell decay spectrum. 

• Large running-coupling effects, which completely dominate the NNLO correction to 
the spectrum [39] if the coupling at NLO is renormalized at rrn,. Large running- 
coupling effects reflect the fact that the typical gluon virtuality is much smaller 
than nib. Naturally, large running-coupling corrections appear also at higher orders. 
Real-emission corrections proportional to Ci?/3g~^as", where /3o (2.16) is the leading 
coefficient of the beta function, were computed in [18] to all orders^. In Sec. 2.3 below 
we will show that the series composed of these terms alone (i.e. with no exponenti- 
ation, no effect of real-virtual cancellation) cannot be considered a viable prediction 

^Another calculation of these corrections is reported in Ref. [29]. There, an infrared cutoff was applied. 
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Figure 1: Comparison between the differential spectra for B — > Xgj decay (through O7 only) 
obtained by DGE, matched to NLO or NNLO, with fixed-order results in space {x space). All 
curves correspond to the same value of the b-quark mass, the PV pole mass nib — 4.89 GeV. The 
end of partonic phase space, = mfj/2, is denoted by an arrow. 



for the spectrum anywhere in the peak region. This series is Borel summable up to 
^ 2 GeV, but ceases to be so above this scale, where the Borel integral diverges 
for u — > 00. 

The numerical results for the spectrum obtained by fixed-order calculations at NLO and 
NNLO, as well as running-coupling corrections beyond this order, are shown in Fig. 1, 
with the pole mass is set^ to = 4.89 GeV. The lack of convergence in the peak region 
is apparent. 

These characteristics of the fixed-order result point towards the necessity for (1) re- 
summation of all large perturbative corrections (2) systematic separation between pertur- 
bative and non-perturbative contributions and appropriate parametrization of the latter. 
Dressed Gluon Exponentiation (DGE) offers a framework to do so, utilizing directly the 
on-shell scheme [17,18,40]. Owing to its inherent infrared safety, the on-shell decay spec- 
trum itself provides then a first approximation to the physical B-meson decay spectrum. 
Fig. 1 shows these results as well. Other approaches to compute the spectrum [16, 27, 29] 
have been based on introducing an infrared cutoff as means of separation between per- 

^Tlie pole-mass value nib ~ 4.89 GeV is obtained by Principal Value Borel summation starting from 
mf^ = 4.20 GeV, as explained in [17]. 
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turbative and non-pcrturbative corrections. With a cutoff in place, a first approximation 
to the B-meson decay spectrum is obtained only upon convoluting the perturbative re- 
sult with a leading-power non-perturbative quark distribution function, or "shape func- 
tion" [25-27,30-33,37,38]. Ref. [17] has shown that a cutoff and a leading-power "shape 
function" can be avoided, and replaced by a resummation of the perturbative expansion, 
a prescription for the renormalon singularities (Principal Value of the Borel sum) and 
power corrections. This involves of course making certain assumptions on the perturbative 
expansion and on non-perurbative corrections. Aside from the separation issue, the calcu- 
lation in Ref. [17] applies a novel approach to resummation per se. Different higher-order 
corrections are considered important in different approaches: 

• Refs. [16,27] emphasize the significance of Sudakov logarithms and utilize a logarith- 
mic accuracy criterion and an infrared cuoff. 

• Ref. [29] emphasizes the significance of running-coupling corrections (dismissing Su- 
dakov logarithms) and utilizes large-/3o resummation with a Wilsonian cutoff. 

• Ref. [17] resums Sudakov logarithms as well as running-coupling corrections in the 
Sudakov exponent and utilizes additional information on the Borel transform of the 
exponent to complement the logarithmic accuracy criterion. 

Upon resumming the perturbative expansion for the matrix elements Gij as a function 
of the cut £^0) normalized by Gij{Einm), Eq. (1.1) can be written as: 

TiB Xsj,E > Eo) = \V,^VCf (mrim,))' ml (1.3) 

X J2 Cf{fi)Cf{ix)Gij{E^i^,^i) 

where it is assumed that .Emm ^ "^6/2 so Gij{E-aan-, IJ) can be computed using fixed-order 
perturbation theory (our choice will be i?min = "),;,/20) while Eq can take any value, in 
particular, experimentally-relevant one: Eq > 1.8 GeV. 

The essential element in the calculation of resummed spectra near threshold, namely 
[Gij (i?o)/Gjj(i?min)] Resummed i^-'^)^ Sudakov factor, which sums up the domi- 

nant corrections owing to multiple soft gluon emission that exponentiate in moment space. 
In DGE the Sudakov exponent is computed as Borel sum, combining Sudakov resummation 
with the resummation of running-coupling effects. This immediately exposes the power- 
like infrared sensitivity of the on-shell matrix elements in the form of renormalons. This 
allows for 

• Explicit cancellation of the leading renormalon ambiguity, the one associated with 
the definition of the pole mass. 

• Definite regularization of all renormalons, using the Principal Value prescription. 

This definition of the perturbative sum amounts to a systematic separation between per- 
turbative and non-perturbative corrections. This procedure uniquely defines the non- 
perturbative Fermi motion effect, distinguishing it from the radiation effect that is common 
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to an on-shell heavy quark and to one that is part of a meson. As already mentioned, 
in contrast with other formulations, no infrared cutoff is needed here. Therefore, Fermi 
motion effects enter exclusively through power corrections. These power corrections are 
power-enhanced at large N: they scale as powers of NA/mi,, so they can modify the per- 
turbative on-shell spectrum by 0{1) corrections near threshold while almost not affecting 
the first few moments and the distribution away from threshold. According to the renor- 
malon structure of the exponent, non-perturbative corrections start at the third power of 
NA/mh, making just a small correction to the spectrum. The numerical significance of 
these power corrections will be analyzed in Sec. 4 below. 

As shown in Fig. 1 the DGE spectrum is qualitatively different from the fixed-order 
results. It is characterized by 

• Approximate physical support: the spectrum smoothly extends beyond the perturba- 
tive endpoint, = mb/2, and tends to zero close to the physical one, = Mb/2. 

• Mass-scheme independence, owing to the explicit cancellation of the pole-mass reror- 
malon ambiguity. 

• Stability in going from order to order, reflecting the fact that the dominant higher- 
order corrections are indeed resummed. 

The DGE B — > Xg^ spectrum, and its first few moments with experimentally relevant 
cuts, were computed in Ref. [17]. Later, when experimental data for the average energy the 
variance appeared [12 14], these predictions where found to be in good agreement [19,41]. 

Recently there has been significant progress in higher-order calculations. The NNLO 
corrections to the spectrum associated with the matrix element Gtj{Eq) of the magnetic 
dipole operator 

07= 3^ mMsi^'^" 5(7^,(1+75)6 (1.4) 

have been computed in full [39,42]. The resummed spectrum of Ref. [17] already included 
NNLL corrections through the Sudakov factor [17,43]; however, it included non-logarithmic 
corrections to NLO only. One of the tasks of the present paper is to match the resummed 
spectrum to full NNLO accuracy in the O7 sector. We also consider other operators in 
the Weak Hamiltonian, which have been so far computed to 0{as) only. Knowing that 
independently of the nature of the short-distance interaction, all important contributions 
in the peak region necessarily involve the same Sudakov factor, we compute resummed 
spectra for individual matrix elements Gij{EQ), determining the hard coefficient functions 
at 0{as) from known results. 

In addition, significant progress towards NNLO calculation of the total BF was re- 
cently made [44,45]: two-loop matrix element of the O7 operator have been computed 
in full. This adds to the already existing NNLO results in the framework of the effective 
Weak Hamiltonian, which includes the matching coefficients at the Weak scale [46], partial 
information on the evolution matrix [47,48], as well as /3oCks^ contributions to B — > X^^ of 
several matrix elements [49]. It is our aim here to set a framework where the state-of-the- 
art calculation of the spectrum can be used together with that of the total BF. This would 
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be particularly important once the NNLO calculation is complete. A detailed knowledge of 
the partial BF as a function of the photon-energy cut E-y > Eq can help making good use 
of the data, since low cuts, such as Eq = 1.8 GeV, are characterized by large systematic 
experimental errors, in contrast with higher cuts, such as Eq = 2.0 GeV, that, in turn, 
requires larger extrapolation. 

An important new ingredient in the calculation of the spectrum that we develop in this 
paper is the use of the analytic structure of the perturbative result in moment space when 
writing the resummation formula. It is a general problem in the application of Sudakov 
resummation, that the region of interest may extend far beyond the asymptotic Sudakov 
regime where the logs are large and the coupling is small. An immediate implication is that 
the rcsummcd result depends on the (often implicit) assumptions made concerning non- 
logarithmic 0{1/N) higher-order corrections. Non- logarithmic corrections arc of course 
included to some fixed order in ct^ in the process of matching the rcsummcd spectrum into 
the fixed-order expansion. This procedure, however, may not be sufficient to avoid bias 
of the result due to the resummation of logarithms in the region where the logarithms are 
not at all dominant. 

The B — > Xg^ spectrum provides an important motivation to address this problem, 
since the region where the logarithms alone dominate is rather small, at least up to the 
NNLO level [39]. To make good use of perturbation theory it is therefore important to 
impose additional constrains on the resummation formula. Such constraints are indeed 
available: the analytic structure of the perturbative result (at any order) in moment space 
is known fairly well: perturbative coefficients are composed of harmonic sums and rational 
functions whose singularities appear on the negative real axis in N space. The rightmost 
singularity appears at = — J where J, a non-negative integer, corresponds to the power 
fall of the spectrum, dV/dx ~ x*^, in the x ^ limit. In B — > AT^T, for most matrix 
elements, J = 3. Thus, there is quite a strong suppression of the spectrum at small x, 
which obviously would not be respected by a generic large-x resummation formula that 
accounts only for ln'(l — x)/{l — x) terms. A power fall is generally expected when the 
a; — ^ behavior is dominated^ by phase space. 

The remainder of this paper is organized as follows: Sec. 2 is devoted entirely to the 
normalized spectrum of the G77 matrix element, corresponding to the magnetic dipole 
operator, O7. G^^ is the only matrix element contributing at C(l), while other Gij start 
at 0{as)- Moreover, G77 is the only matrix element for which the full 0{as^) (NNLO) 
result is available [39,42,44,45]. This facilitates matching of the resummed spectrum to 
NNLO as well as performing an in-depth analysis of the perturbative expansion. In Sec. 2.1 
we briefly summarize the main assumptions and the necessary formulae in the application 
of DGE to the radiative B decay spectrum. In Sec. 2.2 we reformulate the resummation 
formulae under constraints on the analytic structure in moment space, in order to have a 
good description of the small E~^ tail. Next, in Sec. 2.3 we analyze the all-order result in 
the largc-/3o limit. Sec. 3 is devoted to computing the resummed spectra for the matrix 
elements Gij other than G77. In Sec. 4 we compute the total BF and incorporate the 

^It does not apply when this hmit is characterized by singular matrix elements, as is the case of the 
chromomagnetic operator contribution to B — > Xs^. 
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resummed spectra of Sees. 2 and 3 into Eq. (1.3), to obtain the BF as a function of Eq. 
We also study there numerically the relation between renormalons, power correction and 
support properties. Finally, we present predictions for the first few moments as a function 
of the cut and analyze the theoretical uncertainty. In Sec. 5 we shortly summarize our 
conclusions. 

2. Resummed spectrum for the magnetic dipole operator O7 
2.1 Dressed Gluon Exponentiation 

DGE is a general resummation formalism for inclusive distributions that is designed to 
address kinematic threshold problems in QCD where there is interplay between Sudakov 
logarithms, running-coupling effects and parametrically-enhanced power corrections. The 
formalism will not be described here in any detail; we refer the reader to a recent review [40]. 
Here we concentrate on the application to B — > Xs'j [17], briefly summarizing the main 
assumptions and the necessary formulae. 

For simplicity we consider in this section the normalized spectrum associated with the 
magnetic dipole operator only, postponing the calculation of the spectra of other matrix 
elements to Sec. 3, and the calculation of the overall normalization of the BF to Sec. 4, 
where the contributions of the different matrix elements are combined with the appropriate 
Wilson coefficients. 

Formulating Sudakov resummation in moment space 

The normalized differential spectrum associated with the magnetic dipole operator takes 
the form 

dTo^{x) ^ ^ (^asinib)) 6{1 - x) + R {as{mb),x) , 



Tqy dx 



V{as{mb)) = l + CFf2kn f^I^Y 

n=l ^ ' 

R{as{mb),x) = i?3i„g. {as{mb),x) + i?,eg. {as{mb),x) (2.1) 



n=l 

00 / f \\ n 



i?.g. Mmb),x) = Cf f;r-(x) f^f^^l^ 

n=l ^ ' 



where R^ing. contains only plus distributions of the form [ln'(l — x)/(l — x)]_|_ where at 
order ag"', < Z < 2n — 1, while R^^s. is integrable fov x ^ 1. The constants kn are deter- 
mined such that the integral of the normalized spectrum would be exactly unity: 

kn = - f dxr^^'ix), (2.2) 
Jo 

at every order in the perturbative expansion. 
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The partially-integrated matrix element defined with a -E-y > Eq cut, normalized by 
the fully-integrated one, is 



1 dToj{x) 



dx "^U7V-^^ _ ^2.3) 

^77(0,7715) Jx=2Eo/mt ' dx 

Defining the spectral moments of Eq. (2.1) as in Eq. (1.2) one may resum the perturbative 
expansion corresponding to Eq. (2.3) as follows: 



G77{Eo,mh) 



. G77{0,mb) J Resummed 27ri 



1 dN f 

Jc-ioo N -l\ 



2Eo 

nib 



l-N 



+ I dxAR{as{mb),x), 

x=2Eo/mi, 



H {as{mb),N) x S\xd{N,mb) 

(2.4) 



where the first line, written as an inverse-Mellin transform^, is the dominant contribution, 
which contains in particular the contributions originating in V in Eq. (2.1) as well as all 
the log-enhanced contributions to the matrix element, -Rging. in Eq. (2.1), while the second 
line contains some residual real-emission terms that are integrable for x — 1. 

The Sudakov factor 

The Sudakov factor, Sud(iV, m^,), sums up the logarithmically-enhanced corrections origi- 
nating i?sing. to all orders. These corrections exponentiate: 



Sud(A^, rrib) = exp ICp 



TT 



a 



In DGE the exponent is written as a Borel sum: 



Cf 



Sud{N, nib) = exp < / — T{u) 



du 



/3o Jo u 



mt 



Bs{u)V{-2u) 



r(iv) 



(2.5) 



1 



- Bj{u)T{-u) 



T{N) 



1 



T{N - 2u) r(l - 2u) 

(2.6) 



T{N~u) T{l-u) 
where Pq is defined in (2.16) below. Bs{u) and Bj{u) are the Borel transforms 

S{as{fi)) = ^ / duTiu)(-) Bs{u) 

Po Jo / 



J{asin)) = -f / duT{u)(-) Bj{u) 
Po Jo VM / 



(2.7) 



of the Sudakov anomalous dimensions of the quark distribution function [43] and the jet 
function [50], respectively, and T{u) is the Laplace conjugate of the 't Hooft coupling [51]: 

't Hooft ^00 ^ ^2 \ " 



T{u) 



PoOis 



IT 

u5 ^—uS 



.00 /^2\' 



dA 



ln(MVA') = --51n{l + — 



1 



r(i + uS) ' 



A 



din IJ? 



1 



-A\1 + 5A), 



5A 



(2.8) 



''The contour running parallel to the imaginary axis is assumed to be to the right of all the singularities 
of the integrand. 
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with 5 = 

Let us point out that the particular A''-dependence in (2.6) is based on the Melhn 
transform of the plus distributions of the log-enhanced terms, [ln'(l — x)/(l — x)]_^, mak- 
ing no further approximation. Eq. (2.17) in Ref. [17] is based on an approximation to this 
result, a minimal scheme where only \nN terms are exponentiated; it differs from (2.6) 
above by 0{1/N) terms. As discussed in Sec. 3.4 in [52], Eq. (2.6) goes beyond this mini- 
mal scheme in exponentiating a particular class of 0{\/N) terms. In the present paper we 
will take one further step in this direction (see Sec. 2.2 below) and modify the exponent to 
conform with the analytic structure of the full perturbative result in moment space. The 
Sudakov factor then takes the form (2.14). It should be clear that these modification do 
not affect the large-TV behavior, and their influence on the distribution in the peak region 
is small. 

The calculation of the exponent proceeds as in Rcfs. [17, 52]. Formally, Eq. (2.6) 
(or Eq. (2.14) below) is computed with NNLL accuracy, using the known [17,39,43,53- 
55] 0{as^) expansions of 5(as(/i)) and J{aa{lj))^- However, since the Borel integral is 
evaluated, not expanded(!), large subleading corrections, notably running-coupling effects 
are accounted for to all orders [17,52,56-59]. Importantly, the Borel integrand presents 
singularities at integer and half integer values of u. These induce ambiguities that scale 
as integer powers of NK/rriiy and A^A^/m^ that arc inherent to the perturbative quark 
distribution and the jet, respectively. In DGE the perturbative exponent is defined as the 
Principal Value of the integral in Eq. (2.6), while non-perturbative corrections are assumed 
to follow the ambiguity structure of the exponent. 

The details of the spectrum arc dictated by the two functions: Bs{u) and Bj{u) in 
Eq. (2.6). In this paper we adopt the approximations to these two functions that were 
developed and used in previous papers, namely Eq. (3.27) in Ref. [52] (or Eq. (2.35) 
in Ref. [17]). These approximations are based on the known NNLO expansions of the 
anomalous dimensions and on additional constraints on the behavior of these functions 
away from the origin. These are particularly important in the case of the quark distribution 
function Bs{u), since NA/nih is not necessarily small; the contribution of Bj{u) away from 
u = has a rather significant suppression, since Nh? /mf^ ^ 1 at any relevant A^. The 
additional constraints can be briefly summarized as follows [17, 52] (see also the recent 
review [40]): 

• Properties of the large-Po results [18] that are expected to hold in the full theory: the 
Sudakov anomalous dimensions in Eq. (2.7) have no Borel singularities. Moreover, 
their Borel transform vanishes at certain integer values of eliminating some poten- 
tial renormalon singularities in the exponent of Eq. (2.6). Specifically in the quark 
distribution function there is one zero at u = 1: Bs{u = 1) = 0, so 0{N'^Is?' /mf^) am- 
biguities are absent while higher^ power ambiguities are present. This suggests that 
the width of the spectrum is "protected'' from non-perturbative power corrections, 
while higher moments receive such corrections. 

®The two-loop results for the Sudakov anomalous dimensions have recently been checked by additional, 

independent calculations [20,21]. 

®The 0{N K/mb) ambiguity cancels against that of the pole mass [17,18]. 
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• The computed value of Bs{u = 1/2), corresponding to the large-order asymptotic 
behavior of the Sudakov exponent (which is different from its large-/3o hmit). This 
calculation (see Sec. 2.3 and Appendix B in Ref. [17]) is based on the known cancel- 
lation mechanism [18] of the leading, 0{NA/mh), renormalon ambiguity in the expo- 
nent with that of the pole mass, the known structure of this Borel singularity [17,60] 
and the perturbative expansion of the ratio between the pole mass and m^^. 

It was further observed in Refs. [17,52] that given the constraints described above — in 
particular, the expansion of Bs{u) near the origin and its values at n = 1/2 and at n = 1 — 
and so long as Bs{u) does not get large at intermediate values of u (specifically for u ~ 3/2) 
the support properties of the resummed spectrum are close to these of physical spectrum. 
In this scenario power corrections are expected to be small. When making predictions we 
shall not assume that this is necessarily the case, but instead allow for variation of Bs{u) 
and for power corrections as explained below. 

Renormalons and power corrections in the exponent 

We base our analysis on the parametrization of Bs{u) in Ref. [52] (see Eqs. (3.27) to (3.29) 
there), where 

Bs{u = 3/2) = -0.23366 C3/2. (2.9) 

C3/2 = 1 is the default value used in Refs. [17,52], and variation of C3/2 between 0.1 and 10 
is considered. Here, however, we take one step further in the way power corrections are 
taken into account. Introducing power corrections based on the renormalon ambiguities in 
the quark-distribution part of the Sudakov factor (2.14) we have [52]: 



Sud^"^\iV, mb] 



Sud^'^\N,mb) 



X 

PV 



k k=3 

where the A?^-dependence of each renormalon residue is carried by 

r(iv + J) r( J + 1) 



exp|f;er(^)'i?("H^,w|, (2.10) 



R^'^\N,k/2) = Resr(-2ti) 
and where we defined 



T{N + J-2u) r( J + 1 - 2n) 



u=k/2 



(2.11) 



PV ^Cf^ ,pv T{k/2) 



eV^^TTfr Bs{k/2), (2.12) 

where f^^ are dimensionless non-perturbative parameters. Assuming that the power cor- 
rections are of order of the ambiguity itself, Eq. (2.12) with /^^ ~ 1 provides an estimate 
of their magnitude. This assumption will be eventually tested by data. 

The application of (2.10) is limited in practice for several reasons: (1) going to large 
Ej an increasing number of power terms become relevant. It is not clear a priori how many 
would be needed in a given situation; (2) Bs{k/2) with A; > 3 is not known beyond the 
large-/3o limit, and this limit most likely provides an overestimate. If one varies C3/2 in 
(2.9) over a large range, also the magnitude of the power correction €3^ in Eq. (2.12) varies 
a lot if /g^ is kept at its natural size, /J^ « 1. 
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Fortunately, there is another selection criterion that can be used to determine the 
allowed range in the parametrization of Bs{u) as well as in the size of the power correc- 
tions Cfc: the support properties. In the calculation itself the heavy-meson mass is not 
used. However, since it sets the physical support, it can be used to distinguish acceptable 
spectra from non-acceptable ones. In Sec. 4.3 we perform a numerical analysis examining 
the parameter space of C3/2 and f^^ under constrains on the support of the corresponding 
spectra, see Eq. (4.12) and Fig. 9 there. 

Matching the resummed spectrum to fixed order 

Given the Sudakov factor, one can determine the hard coefficient function H {as{mb),N) 
as well as the residual terms AR{as{mi,) , x) in Eq. (2.4) order-by-order in perturbation 
theory from the expansion of the differential spectrum in Eq. (2.1). It should be noted 
that as far as the terms that are regular for x ^ 1 are concerned, the separation between 
the contributions that are taken into account in moment space and the residual terms 
that are included directly in x space is arbitrary. Moreover, the matching between the 
resummed exponent and the fixed-order expansion can be done in a variety of ways that 
differ by subleading corrections. In Appendix A we derive explicit expressions to ©(a^^) 
for H (04.(717.5), N) and AR{as{mi,), x) based on the available NNLO results of Rcf. [39,42]. 
In doing so we rely on previous experience in the applications of Sudakov resummation to 
inclusive distributions in QCD, e.g. [17,52,56,57,59], in the following ways: 

• we give preference to moment space, so AR{as{m,h), x) reduces to small corrections 
that vanish at a; ^ 1. This guarantees, in particular, smooth transition from the 
perturbative region x < 1 to the one above x = 1. 

• in moment space we use "log-R" matching [61], where the perturbative expansion of 
the logarithm of the spectral moments is constructed as a sum of the perturbative 
expansions of In Sud(iV, m(,) and of InH {as{mb), N). Consequently H {as{mh), N) 
itself is constructed as an exponential function. 

Following these considerations we arrive at the NNLO matching formula of Eq. (A. 28). 

2.2 Analytic structure in moment space and the small— Ej asymptotics 

In general, Sudakov resummation is a parametrically-controlled approximation in a specific 
kinematic region where the logarithms are large and the coupling is small. The logarithmic- 
accuracy criterion applies only if the product of the coupling times the logarithm is small 
enough. The application to b decay, where the coupling at the hard scale is ct^ (771-5) — 
0.2, is a borderline case at the outset: near threshold, where the logarithms are really 
large, the logarithmic-accuracy criterion does not hold and non-perturbative corrections 
become important, while away from threshold the logarithms are no more the dominant 
perturbative corrections. 

As discussed above, DGE is designed to deal with the first problem, extending the 
applicability of the resummation deeper into the threshold region, x ~ 1. In this section 
we would like to address the second, namely the application of the resummed spectrum in 
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the region where the logarithms are not necessarily dominant. Usually, Sudakov-resummed 
spectra suffer from artifacts when evaluated away from the Sudakov region. We will show 
that by imposing constraints on the analytic structure of the Sudakov factor in moment 
space, one can extend the applicability of the resummed spectrum away from the Sudakov 
region, and even use it for x ^ 1, where in the absence of such constraints artifacts from 
the resummed ln^(l — x)/{l — x) terms are significant. 

Much like the problem, the solution we suggest is general. We nevertheless consider 
here the concrete case of the photon-energy spectrum in b ^ sj. As we shall see, it 
provides a particularly good example because of the strong suppression of the differential 
spectrum at small x, which makes potential resummation-artifacts more pronounced. 

As was observed by Melnikov and Mitov [39], at small x, the differential 6 — 57 
spectrum dTo-,/dx falls as x^: 



This is a general property that holds to all orders in the pcrturbative expansion. In 
particular, the known two-loop results [39,42] as well as all-order results in the large-/3o 
limit [18] — see Appendix B.3 below — all share this cubic suppression. This suppression is 
partially owing to the fact that the available phase space shrinks as E^dE^ for E^ — 0, and 
partially a consequence of the dynamics: the coupling of the photon to the flavor-changing 
current in (1.4) involves F^^ = d^A^ — d'^A^, which translates into a power of the photon 
momentum compared to the usual coupling. In the squared matrix element this results 
in two powers of E^ on top of the phase space suppression, namely dVoj ~ E^dE^. 

It is obvious that any resummation procedure that considers only the terms that are 
singular at x —>■ 1, namely, Q;s"ln'^(l — x)/{l — x), is bound to generate artifacts in the 
region where these terms no longer dominate. The largest artifacts are associated with 
the most subleading logarithms, 1/(1 — x), which behave as a constant for x — > 0. Such 
resummation artifacts arc limited in size by virtue of matching the resummed Sudakov 
factor to the fixed-order expansion. However, matching to a fixed order may not always be 
sufficient to guarantee a good approximation away from the Sudakov region, especially in 
cases where the coupling is large, so subleading perturbative corrections are not negligible. 

The strong suppression of the differential spectrum in Eq. (2.13) makes such artifacts 
particularly important: if only a,5"ln^(l — x)/{l — x) terms are resummed, resummation 
artifacts that behave as a constant at small x are expected to appear at any order. NNLO 
matching guarantees that such terms only appear at 0{as'^) and beyond. Nevertheless, at 
sufficiently small x these can easily compete with the terms in Eq. (2.13), and eventually 
dominate. Consequently, the resummed spectrum would develop a (small, 0{as^)) constant 
X ^ tail instead of the correct cubic fall-off. Obviously, this should be avoided. In the 
following we develop a formalism where such artifacts are systematically avoided to any 
order by imposing the correct analytic structure on the moment space Sudakov factor. 

The basic observation underlying our method is that there is a direct correspondence 
between the small-x asymptotic behavior and the analytic structure in moment space 




(2.13) 
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Figure 2: Comparison between the differential spectra for B — > Xgj decay (tlirough Oj only) 
obtained by DGE using the exponent (2.14), matched to NNLO according to (A. 43) with J = 
through 4 (and fi = rrib), and results obtained in x space with NNLO using Borel summation of 
running-coupling effects in the SDG approximation according to Eq. (2.23). In the plot on the 
left the horizontal axis ends at the physical endpoint Ej = Mb /2 while the partonic endpoint is 
denoted by an arrow. The plot on the right hand side enlarges the tail of the distribution; it ends 
at the partonic end of phase space. 



(where moments are defined according to Eq. (1.2)): a pole at = — J corresponds to a 
small-x fall-off of the form: dV/dx ~ x"^ . 

Motivated by Eq. (2.13), we assume the asymptotic behavior dT/dx ~ x"^ (in om' case 
J = 3) and construct the Sudakov-resummed spectrum to accommodate this behavior. In 
order to capture the x'^ suppression at x ^ at the level of the Sudakov exponent we 
modify Eq. (2.6) into: 



Sud (A, mfe) = exp <^ — / — T[u) 

Po Jo ^ 



ml 



Bs{u)V{-2u) X (2.14) 



r(iv + j-2ti) r(j + i-2n)y ^ ^ >\y{n^j-u) r(j + i-u) 

This guarantees that no poles are generated for N > — J and therefore already before any 
matching is done the Sudakov factor would not give rise to a tail that falls slower than x'^ . 

In Appendix A. 3 we develop the matching procedure of Eq. (2.14) with the NNLO 
results, in analogy with what was done in Appendices A.l and A. 2 for the J = case. 
Similarly to Eq. (2.14), the matching coefficients are constructed under a constraint on the 
analytic structure in moment space: no poles should appear for N > — J, and so the small- 
x asymptotic behavior would coincide with that of the fixed-order result, dT/dx ~ x'^ . 
The final matching formula, for a general J, appears in Eq. (A. 43). Note that while 
both the exponent (2.14) and the matching coefficients entering (A. 43) vary with J, the 
NNLO expansion does not; only higher orders do. Moreover, the log-enhanced terms are 
J independent, only 0{1/N) terms depend on J. 
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Figure 3: The partial BF ratio, BF^E^ > Eq)/BF{Ej > 0), corresponding to magnetic dipole 
operator, as a function of the cut Eq, computed by DGE (2.14) with J = 3, matched to NLO 
(dashed) and NNLO (full line) according to Eq. (A. 43). In both cases the renormalization scale 
in the matching coefhcient is fi = nib. The shaded area represents the uncertainty band at NNLO 
while the two external dashed lines represent the NLO one. The lower NLO uncertainty band trails 
the central line for the NNLO matching prediction for Eq < 2.15 GeV. The width of the uncertainty 
band is reduced roughly by a factor of two by including the NNLO corrections. The calculation of 
the theoretical errors is explained in Sec. 4.2 below. 



The effect of varying J is shown in Fig. 2. As expected, variations in the peak region as 
a function of J are small, while variations in the tail are rather significant. In particular, as 
shown in the plot on the right, the J = curve approaches a constant at small x, the J = 1 
one falls as x^, etc. These are artifacts of the resummation. Such artifacts are absent for 
J = 3 (or larger) where the power fall-off is determined by the perturbative expansion itself, 
Eq. (2.13). Obviously, the natural choice in our case is J = 3: singularities at = —3 do 
appear in the perturbative expansion, so they should be allowed in the resummed spectrum. 
We will use the J = 3 spectrum as the default choice below. 

Having set the final resummation and matching formulae, let us examine the predic- 
tion for the partial BF as a function of the cut. The result is shown in Fig. 3. The figure 
compares the results obtained at NLO and at NNLO. Note that the Sudakov factor, which 
is computed as a Principal Value Borel sum, is the same (formally it has NNLL accuracy) 
— only the matching coefficient (A. 43) is truncated at different orders. The error is sig- 
nificantly reduced at small Eq cuts, where the corrections in the matching coefficient are 
essential. 
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2.3 All order results in the large /3o limit 

Running-coupling effects often constitute a major part of the radiative corrections in 
QCD [62-64]. The typical situation is that the average gluon virtuality is lower, some- 
times quite significantly lower, than the hard scale in the process, m. Thus, when the 
latter is used as the default renormalization scale for coupling {fj, = m), one encounters 
large corrections involving powers of {(3oOis{m)/TT), to any order in perturbation theory. In 
simple cases perturbative predictions can therefore be significantly improved by using the 
BLM scale instead. Differential spectra involve other parametrically-large corrections, such 
as Sudakov logarithms. Since they inherently depend on several scales, a straightforward 
scale-setting procedure is excluded. In such cases resummation is necessary. 

In case of the b Xgj spectrum 0{CFf3oas'^) running-coupling (BLM) corrections 
have been computed long ago [65]. When using /i = ruh, these corrections are very sig- 
nificant. A couple of years ago running-coupling terms in the G77 spectrum, CFPo~^Ois"' , 
were computed to all orders [18,29]. Recently, the full ©(a^^) calculation of G77 was per- 
formed [39,42], finding that the additional, non-BLM corrections are, instead, moderate. 
Figs. 4 and 1 present these results for the normalized differential 6 Xgj spectrum. 
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Figure 4: Convergence of the perturbative expansion in Eq. (2.19), generated by running-coupling 
effects only, for the O7 decay spectrum: partial sums of order ©(a^") with n = 1 through 4 are 
shown together with the Borel sum of Eq. (2.15). Left: the whole range of energies on a logarithmic 
vertical axis. Right: a restricted range on a linear scale. In both plots the horizontal axis ends at 
— mf,/2, as dictated by the support of these functions. The d{E^^ — mf,/2) is not shown. In the 
plot on the right hand side the additional full lines at NNLO and beyond include the full 0{as^) 
corrections based on the calculation of Ref. [39], see Eq. (2.22) and (2.23) below. 
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The plots clearly show that the Cf(3q~ cks" corrections are very large. In particular, if 
the expansion is performed in terms of ^^(m;,), the ©(Cir/^octs^) term completely dominate 
the 0{as^) correction. This is true even at fairly large E^, despite the presence of formally 
leading Sudakov logarithms with other color factors. With this numerical evidence it is easy 
to reach the conclusion that Sudakov resummation is not needed, or cannot significantly 
improve the description of the spectrum (see e.g. the discussion in Ref. [29]). The purpose 
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of this section is to demonstrate that this conclusion is wrong, and that multiple soft and 
collinear gluon emission does in fact play an important role in shaping the spectrum through 
higher-order corrections. To this end we compare the fixed-order results with (1) DGE, 
where Sudakov logarithms are taken into account through exponentiation in moment space; 
and (2) the single-dressed-gluon (SDG) approximation, which accounts for real-emission 
running-coupling CF0l^~^as^ corrections, to all orders, but neglects Sudakov logarithms. 
We will show that while the SDG series is Borel-summablc for E-y ^ 2 GeV, it docs not 
provide a viable description of the spectrum anywhere in the peak region, in sharp contrast 
with the DGE result. 



The single— dressed— gluon approximation 



The all-order calculation of the G77 decay spectrum in the large-/3o limit was performed 
in Ref. [18]. Using the scheme-invariant Borel transform with A defined in the MS scheme 
the result takes the form: 



1 dV, 
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To-! '^^ 
with mi, the bottom pole mass, 

and 



6(1 -x) + 



Cf 



duT{u) 



A 



mt 



2 \ " 



B{x, u) 



(2.15) 



+ 



/3o 



11 
12' 



Nr. - ^Nf 



B(x,u) = eh^^^x^l-xy da a (1-a)-" 



(2.16) 



(2.17) 



(1 - xay 



l-4a + a^ 



5 smnu n . 
es" x-^{l -x) " 



(!-«)- 
1-u 

1 



+ 



a 



+ 



+ 



(1 — a;a)(l — u) 1 — a; 1 — xa ^ — x 
1 



+ 



TTU 



(1 - 4x + X^) 



l-a;(l-n)(2-'u) 

2{l-x)^ ' 



1 1 

+ 



1—x 1—u 



+ 



(1 — u)'^x^ 



{l-4x + x^) 1 (x + l)(x2-3x + l 



2Fi([l,l],[2-u],x) 
2-x) 2(1 -x) 



x^ 1 — x (1 — u)x^(l — x) (2 — ■u)x2 (1 — ■u)^x^ j 

In Ref. [18] the Borel function was expressed as an integral over a single Feynman parameter 
(the first expression above) . Here we performed this integral and wrote the Borel function 
in terms of a hyper geometric function. 

Strictly within the large-/3o limit T{u) = 1. It is straightforward to take into account 
running-coupling effects beyond this level by an appropriate choice of the function T{u) in 
the Borel integral^. To account for Pi terms to all orders it is sufficient to use the T{u) of 



^Of course, this can be done in a renormalization-scheme invariant way only as far as two-loop effects 
(/3i) are concerned. 
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Eq. (2.8), see Ref. [51]. This will be done in the perturbative expansion, Eq. (2.19), and 
in the numerical analysis that follows. 

Let us now examine the result of Eq. (2.17). The first observation is that this Borel 
function has no renormalon singularities. To see this note first that (despite appearance) 
there is no double pole at w = 1 inside the curly brackets: for u — > 1 the hyper geometric 
function reduces to 2-^1 (^[Ij 1]) [2 — u],x^ — > 1/(1 — x), so the double pole cancels. Now, 
since the curly brackets contain just simple poles at integer values of u, upon taking the 
svnTiu factor into account, one concludes that there in no renormalon in B[x,u). 

The absence of Borel singularities may be interpreted as an indication that the per- 
turbative expansion in x space should be good. As we shall see this is not the case. First 
of all, at large x the effective scales are paramctrically lower than [18]: the dynamics 
is dominated by momenta of order mj-y/l — x, the jet mass, and mij(l — x), the soft scale 
associated with the momentum distribution of the b quark. The analytic expression in 
(2.17) immediately indicates the presence of the former through its (1 — x)~" dependence; 
to isolate the latter it is convenient to first apply the identity: 



2Fi([l,l],[2-u],x) ={l-u) --2Fi([l,l],[l + u],l-x)+^^{l-x) 



-u^u-1 



(2.18) 

in (2.17), revealing the dependence on (1 — x)~^". Moreover, the absence of renormalons 
does not imply the existence of the Borel sum: the Borel integral (2.15) may not converge 
at u — 00. It turns out that the Borel integral exists for sufficiently small x values, but it 
does not exist at large x. We shall return to discuss this issue below, after considering the 
perturbative expansion order by order. 

As usual, renormalons do appear upon taking moments [18]: when integrating over all 
values of x, including the endpoint region x ^ 1, soft gluons necessarily contribute and 
consequently, infrared renormalon singularities at integer and half integer values of u are 
generated. We emphasize that this is not a disadvantage of the moment space formulation, 
but rather its strength: renormalon singularities are a useful tool to understand the role 
of non-perturbative power corrections. 

Perturbative expansion in x space in the large— /3o limit 

The perturbative expansion of Eq. (2.15) is: 
1 dTor . . Cf °" 



where 
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where the dots stand for higher-order terms, ©(a^"^^) and so on, containing powers of 
Pi/(3q. Up to 0{as^) such terms do not appear, and Eq. (2.19) reduces to 
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TT 
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The coefficients can be obtained order by order from the expansion of ^B(x,u) = 
X^^o "i^n^ix) u"'/n\ in Eq. (2.17). To this end we need the expansion of the hyper geometric 
function, which is available based on Ref . [66] . We present this expansion and the resulting 
coefficients up to 0{as'^) in Appendix B. 

It is straightforward to include the remaining Cp and CaCf terms at 0{as^) based 
on the results of Ref. [39]^: 
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The partial sums of Eq. (2.19) as well as the Borel sum of Eq. (2.15) include the large-/?o, 
(x) terms at each order. Thus, to have a complete result to 0{as'^) we now include 
the non-BLM {Cf and Ca) terms in the square brackets of Eq. (2.22). When using Borel 
summation for the running— coupling contributions the normalized spectrum at NNLO can 
be expressed as: 
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The results are shown together with the pure running-coupling contributions in the plot 
on the right hand side in Fig. 4. As already concluded in Ref. [39], these additional 
contributions are moderate, so at least at this order, the running-coupling contributions 
dominate. 

Let us now examine the convergence of the expansion and the possibility to sum up 
the series a la Borel. Fig. 4 summarizes the numerical results using ruf, = 4.875 GeV and 
A = 0.332 GcV with Nj = 4 and two-loop running coupling. The partial sums obtained 
by truncation of Eq. (2.19) at increasing orders converge well at small energies ^ 1.5 
GcV and less so as the energy increases. Nevertheless, the increasing-order contributions 
do decrease monotonically. This decrease is related of course to the absence of renormalons 
in Eq. (2.17). This usually means that the series is Borel summable. However, since it is 
a divergent series there is no straightforward relation between partial sums and the Borel 
sum. 

As anticipated, the series is in fact Borel summable only below a certain energy, 
X < Xraax.'- the Convergence of the Borel integral at large u is guaranteed only owing to 



®We wish to thank Melnikov and Mitov for providing us with their final expressions in the form of a 
Maple file. 
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the suppression by (A^/m^)**; as already observed in Ref. [18], because of the presence 
of (1 — a;)~^" contributions in B{x,u) (see Eq. (2.18)), it is predominantly the soft scale 
mi){l — x) which sets the argument of the couphng at large x. Based on these considerations 
we can deduce from Eq. (2.17) an estimate of where the Borel sum will cease to exist: 

^max ^'^(^_ e5/6^/^^j ^ 2.06 GeV. (2.24) 

Indeed, as shown in Fig. 4 the Borel sum is close to the high-order partial sums 
at small energies and to lower-order ones at higher energies^, until it eventually reaches 
a peak, bends downwards and becomes negative. Soon after it breaks down completely 
owing to non-convergence of the u integral at u — > oo, in accordance with Eq. (2.24). 

We note that increasing-order partial sums become quite different from each other at 
Ej ^2 GeV; this obviously means they cannot be thought of as an approximation to the 
spectrum. Also the corresponding Borel sum, although unique where it exists, does not 
look anything like a physical spectrum would. Its complete breakdown according to (2.24) 
is indicative of the low scales involved. But even if the bottom mass were much higher, this 
SDG approximation cannot be expected to describe the Sudakov region. Running-coupling 
corrections are not the most important corrections for Ej — > mb/2. In the large-/3o limit 
the most singular terms at large x are Ug"' Cp Po^~^ — x)/{l — x), while the full 
perturbative expansion contains double logarithms of form a^"" Cf"" ln^"'~^(l — x)/(l — x), 
as well as other subleading logarithms which are associated with multiple soft and coUinear 
emission. With increasing order, the latter become more important at large x compared 
to the former. Finally, returning to Fig. 1, we observe the remarkable difference between 
increasing-order partial sums and the DGE result in the peak region. The gap between 
them should be understood as the contribution of multiple soft and collinear gluon emission, 
with a particular regularization of the divergent sum. In contrast with the fixed-order 
partial sums and the SDG Borel sum, the DGE result provides a useful approximation to 
the spectrum in the peak region. 

3. Resummed spectra for individual matrix elements other than G^j 

In Sec. 2 we considered in some detail the calculation of the normalized spectrum of the 
magnetic dipole operator. It is well known that in contrast with the total width, the 
spectrum is not so sensitive to the details of the short-distance interaction, and therefore 
the G77 spectrum computed above can be thought of as an approximation to the physical 
spectrum. However, to make precision estimates of the partial BF as a function of the cut, 
it is important to compute the spectra of other matrix elements as well. Since fixed-order 
results for all matrix elements but G77 are available to 0{as) only, the accuracy we can 
hope to achieve in the description of the spectra of individual Gij matrix elements is not 
as high. 

^For very low Ej <, 0.3 GeV the Borel sum is still larger than the 0(as'') partial sum; as the energy 
increases it gradually approaches the 0{as^) partial sum, crossing it around E-y ~ 1.65 GeV; around 
Ej ~ 1.9 GeV it crosses the 0{as'^) partial sum and soon after turns over and starts decreasing; finally 
close to 2 GeV it crosses the leading-order result and continues to decrease. 
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We will make the approximation where the matrix elements Gij (for any i and j 
other than 77) are written as 0{as) hard coefficient functions, which are 0{N^) for Gi^ 
interference terms and 0{N~^) for other terms, times the same Sudakov factor discussed 
in the context of G77. This approximation is motivated by the following observations: 

• Independently of the nature of the short-distance interaction, all non-integrable real- 
emission corrections at any order, namely corrections to the moments that scale as 
iV° times logarithms at large AT, are necessarily associated with soft and coUinear 
radiation around the same Born-level configuration involving a b quark in the initial 
state and an unresolved quark jet in the final state. 

• All operators mix into O7 under renormalization. Moreover, the 0{as) result of Gij 

is dominated by the virtual diagrams, which are proportional to the tree level G77. 
This observation was already used in Ref. [49] arguing that 0{Nfa'^) corrections can 
be well approximated by computing both real and virtual diagrams for O7 and only 
virtual diagrams for Oi,2,8- 

It should be noted that while the 0{N^) interference terms do indeed involve the same 
Sudakov factor, 0{N~^) contributions associated with integrable bremsstrahlung correc- 
tions, do not share this property, and in general involve a new jet function. In this respect, 
our calculation of the resummed spectra of Gij that do not involve O7 is a cruder approxi- 
mation; as we shall see the relative weight of these contributions, especially in the Sudakov 
region, is small. 

The structure of this section is as follows: in Sec. 3.1 we discuss the small E"^ behavior 
of the different matrix elements and then extend the matching procedure of Sec. 2 accord- 
ingly, using the known 0{as) coefficients. In Sec. 3.2 we present numerical results for the 
resummed spectra of individual matrix elements, based on this 0{as) matching procedure. 

3.1 The small asymptotic behavior and matching at 0{as) 



Looking at the small photon-energy limit, 
behavior: 



0, we find the following asymptotic 



dGij{Eo,mb) 
dEo 



except for 



= o{Ely, 



(3.1) 



As discussed in Sec. 2.2, the phase-space suppression for £^.y — is E^dE^, and any 
additional suppression or enhancement of the Gij spectrum in this limit depends on the 
dynamics, and therefore on the operators Oi and Oj. To understand the difference between 
Og and O7 in this respect, note that in the Og amplitude the gluon couples to the flavor- 
changing current through G^'^ while the photon is emitted as bremsstrahlung, whereas 
in the O7 case it is the other way around. While the O7 coupling (1.4) gives rise to a 
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linear dynamical suppression with the photon energy, photon bremsstrahlung involves a 
soft singularity, namely a enhancement. Combining this behavior of the amplitudes 

with the phase-space factor one immediately obtains the different Ej ^ limits in (3.1) for 
G77, Gs8 and the interference Gyg. At 0{as), the behavior of matrix elements containing 
Oi,2 is similar to that of O7. The reason for this is that the photon cannot be emitted as 
bremsstrahlung, but must instead couple to the virtual charm propagating in the loop (the 
diagram with just gluon coupling to the charm loop vanishes). 

Thus in most cases there is a cubic power suppression of the spectrum for extremely 
soft photons, — 0. When matching Sudakov resummation for the high photon energy 
endpoint {Ej = mf,/2) into the fixed-order results, it is therefore useful to take J = 3 as for 
the O7 contribution. As explained in detail in Sec. 2.2, in this way artifacts of resummation 
away from the hard photon endpoint are avoided (the more pronounced soft-energy tail of 
Og and its interference with O7 will be accounted for at fixed order). 

In general, the matching involves a separation of the real-emission contributions be- 
tween momentum and moment space (see Sec. 3 and Appendix C in Ref. [17]). The 
separation we use is defined such that the leading contributions at small A are always 
treated in moment space: 

4-(A) = 0(A2). 



(t>ij (A) = fiij (A) + (A) with 
Introducing the Mellin transform of fiij{A) as in [17]: 



jlij(N) = dxx 



dA 



A=l-x 



pc+too 
J c—ioo 



dN (1 - A)i 



-N 



we get: 



Gij{EuAn,fnb) = Cp 



Ois{mb) 



TT 



^§j=7Re { Ti + In — 



+ C^^£i!I^4.(A) + 0(0 



-L 



TT 



c+ioo 



dN (1 - A) 



l-N 



2'Ki N-1 



TT 



2'Ki N-1 



+ %(A) 



X Sud^'^\iV, mb) 



(3.2) 
(3.3) 



TT 



(3.4) 



where in the second expression used Eq. (3.3), included the Sudakov factor resumming 
large logarithms at ©(a^^) and beyond, and defined 



y«(iV)^^<^,=7Re{r. + 7l?^ln!^}+A., 



(N). 



(3.5) 



Note that fLij{N) are all 0{1/N) (Appendix C in [17]) so v[j\n) is of 0(1) at large for 
j = 7 (the interference terms for any i are 0{1) at large N) while it is 0{1/N) otherwise. 

For 77, 78 and 88 matrix elements the NLO contribution is treated entirely in moment 
space: %(A) = 0ij(A) and iij{A) = 0, so 



{N 



(3.6) 
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and 

~ /An ^2121 11 



9A^ 9iV-l 9iV + l 



1..^4*(iV)-i7.) (3.7) 



11 112 1 11 1 1 



18N + 2 6N 9N-136N + 1 m'^ 18(l + iV)2' 

Note that the latter expression has a pole at AT = 1, owing to the soft-photon singularity. 
However, this should not affect the calculation of Gij{Eiam, mi,) with non-zero Emm so long 
as the contour in Eq. (3.4) is to the right of iV = 1. 

For all the other contributions, namely the ones associated with the operators Oi and 
O2 and their interference with O7 and Og (this includes (^22, 027, ^ii, 4'12 <Pn, 0i8 and 028) 
we introduce a separation according to Eq. (3.2) with a non-trivial |ij(A). The specific 
separation used in Appendix C in [17] was based on defining ??ij (A) as the leading-order 
term in the expansion of 0ij (A) at small A. This simple procedure is not adopted here, 
since it would introduce a constant differential spectrum at very small E-y^ in contradiction 
with the power suppressed behavior of Eq. (3.1). Instead wc will define fjij{A) such that it 
would capture both the small A limit and certain features of the A ^ 1 limit. As shown 
in Appendix C in [17] for small A, 4>ij{A) is C(A), so according to Eq. (3.2) wc require 

for small A: 77ij-(A) ~ (/)^(A = 0) x A + ^(A^). (3.8) 

As explained above (see Eq. (3.1)) for A — 1 the differential spectrum falls as the third 
power 1 — A so the integrated spectrum approaches a constant as the fourth power of 
(1 — A). We therefore require: 

for large A : fjij (A) ~ 0^,- (A = 1) + C» ((1 - A)-^+i) , (3.9) 

where J = 3, such that (1 — A)^ contributions appear through both the moment-space 
expression (fjij) and the momentum-space residual (^ij), while ones with (1 — A)'* for 
h = 1,2,3, which are absent in the physical spectrum, do not appear in either. Both 
^^^•(A = 0) and 0ij(A = 1) can be exphcitly computed from the expressions in Eq. (C.l) 
in [17]. To accommodate (3.8) and (3.9) we define: 



%-(A) = 0ii(A = l) 



1-(1-A)^+^ (i + (J+1)a) 



+ 0^(A = O)(l-A)^+iA. (3.10) 



In computing Gij{Eo,mh)/Gij{Einm,'m,i,) using Eq. (3.4), ^ij(A) is readily obtained using 
^ij(A) = 4>ij{A) — ?7,j(A), while jj-ij{N) is computed according to Eq. (3.3) by the Mellin 
conjugate of the derivative of Eq. (3.10), namely: 

(N) = r dA (1 - A)^+^-' Uj{A = 1) (J + 1)(J + 2)A + <p',JA = 0) (l - ( J + 2)A) 
Jo 

(J + l)(J + 2) N-1 

■bij{A = 1)+ , T\rAT , T I = 0). (3.11) 



'0 



{N + J){N + J+1)^''' ' (Ar + J)(iV + J + 1) 

As required, the first moment ilij{N = 1), corresponding to the total width, is given 
by (pijiA = 1), while the large- A*" asymptotic behavior is determined by (p[j{A = 0). 
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Figure 5: Left: the contribution to the spectrum from the separate operators and their mixing, 
together with their sum, the total width, based on Eq. 4.8. We notice that the dominant contribu- 
tions (from 77 and 27) have similar shapes in the Sudakov region, although their tails for small Ej 
are significantly different. Right: the spectrum of each matrix element Tij (see Eq. (4.8)) has been 
normalised to better display the various shapes. Since effectively Oi = g02j the curves for 17 and 
27 are identical. The same is true for the groups (18, 28) and (11, 12, 22). 



Importantly, when using the matching procedure of Eq. (3.4), where Sud (A^, m^) is given 
by Eq. (2.14) and jlij{N) by Eq. (3.11), we are guaranteed that no spurious singularities 
at A'^ > — J would appear (singularities do appear for N = — J, N = — J — 1 etc.), which 
could alter the small-S^ behavior. 

3.2 Resummed spectra for individual matrix elements 

In the previous sections we established a procedure for matching the Sudakov-resummed 
spectrum with the known fixed-order expansion. For G77 this is done to NNLO, 0(0^^), 
while for the matrix elements of other operators and their interference with O7, to 0{as)- 
This allows us to study the shape of each contribution and eventually have better control 
on the cut dependence of the total BF. 

Figure 5 shows the contributions to the differential BF from each matrix element 
Gij according to Eq. 4.8 below, as well as a comparison between the shapes (normalized 
spectra) of the different matrix elements. In general, it confirms the lore that the shapes of 
the spectra of the most contributing operators do not vary much, and the shape of the total 
spectrum is roughly the same as the one for O7. However, the shapes are not completely 
identical and the details do depend on the relative contributions of the different operators. 
These depend of course also on the coefficient functions. Numerically, the most important 
effect in the experimentally accessible energy range is that of G27 which is comparable 
in magnitude to G77, but is somewhat narrower. The effect it has on the differential 
distribution along the tail, e.g. between 1.8 and 2 GeV, is significant. The G78 spectrum, 
on the other hand, is somewhat wider and it has a higher small-x tail. A rather different, 
and significantly wider shape is presented by G28 and even wider still by G22 (and likewise 
those corresponding to Oi instead of 02)- These affect the differential spectrum at smaller 



-24- 



energies, below 1.8 GeV. 



4. Partially— integrated BF and moments with a photon— energy cut 

In the previous sections we computed the B — > Xg^ spectra using DGE. We improved the 
calculation in Ref. [17] in various ways: we included the NNLO corrections in full in the 
Gji sector (Sec. 2.1 and Appendix A), introduced a new method to constrain the behavior 
of the resummed spectrum away from the Sudakov region (Sec. 2.2 and Appendix A. 3), 
and matched the resummed spectrum with the fixed-order expansion, separately for each 
matrix element Gij (Sec. 3). In this section we will make use of these advancements to 
make predictions for measurable quantities: the partial BF and the first few moments as 
a function of the cut on the photon energy. Wc begin in Sec. 4.1 by computing the total 
BF and discussing its theoretical uncertainty, then in Sec. 4.2, we consider the partial 
BF with a cut and discuss cut-dependent uncertainties, and finally in Sec. 4.3 we present 
results for the first few moments and analyze numerically the relation between renormalon 
contributions, power corrections and support properties. 

4.1 The total BF 

The calculation of the total BF presents a challenge in its own [7,8]. First, this calculation 
has to accommodate a large hierarchy of scales M]Y,mt ^ rrih. This is dealt with [4-8] by 
the well-established machinery of the effective Weak Hamiltonian, which facilitates the re- 
summation of large logarithms, Inmw/frib, to all orders in the strong coupling. Estimating 
the BF therefore involves the perturbative calculation of the matching coefficients at the 
high scales, evolution of the different operators from Mw to nif, and finally, the contribu- 
tion of the matrix elements of the different operators to B — ^ Xsj. The main challenge 
here is at the perturbative level: even the matrix elements present little sensitivity^^ to 
non-perturbative effects. 

The NLO calculation was completed 5 years ago (see [8] and Refs. therein) and some of 
the ingredients have been recently brought to the NNLO level. This includes, in particular, 
the matching coefficients at the Weak scale [46], partial calculation of the evolution ma- 
trix [47,48] and, most importantly, the two-loop matrix element of the O7 operator [44,45]. 
However, other essential ingredients are still missing. Amongst these are the matrix ele- 
ments of other operators, especially G27, which contributes to the BF almost as much as 
G77, despite the fact that its perturbative expansion starts at 0{as), rather than at 0(1). 
The renormalization scale of Gij other than G77 is obviously not fixed at this level. 

Beyond G77, which is known to NNLO in full [44,45], the virtual Pocts^ corrections to 
Gij have been computed in Ref. [49]. It is likely that these running— coupling corrections 
constitute a major part of the NNLO contribution — this has been supported by the finding 
of Ref. [44]. Here we shall make use of this available, partial NNLO information on the 
matrix elements, which are the most important ingredient, neglecting NNLO effects in the 
Wilson coefficients at the Mw scale and in the evolution. A more complete NNLO BF 

possible exception is the large distance sensitivity associated with cc pairs. 
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estimate is due when some of the missing ingredients — in particular G27 at ©(a^^) — 
become available. 

Another difficulty in the evaluation of the BF is its sensitivity to the bottom mass, 
which is raised to the fifth power in (1.1). Amongst these five powers, two are associated 
with the operator itself, and therefore correspond to the short-distance mass m^^, whereas 
the additional three result from phase-space integration, and therefore correspond to the 
pole mass, mfe. As discussed in the introduction, the pole mass has a leading renormalon 
ambiguity at u = 1/2, implying that the perturbative expansion of the total width in the 
on-shell scheme is badly divergent. Writing (1.1) order by order, 

T{B X,j,E> Eo) = |Ftb^t:i' (mrirn,))" ml (4.1) 



327r4 



F 

the expansion coefficients fnilj) are therefore expected to grow factorially with the order n 
owing to running-coupling corrections. Judging from what we know of the pole mass — see 
Appendix B in Rcf. [17] — this divergence sets in already at the first few orders. A similar 
situation occurs for other inclusive decays, for example for the charmless semileptonic decay 
B — >XJu [34-36]: 



"'total {B > Xulv) 



1927r3 



1 + si \-s2W H 

TT \ TT 



(4.2) 



In this case the coefficients have been computed to NNLO in Ref . [70] , where large running- 
coupling corrections already appear. Moreover, in this case it is possible to interpolate 
reliably [52] (sec also Ref. [71-73]) between the first few orders and the known asymptotic 
behavior of the scries, which is set by the n = 1/2 renormalon of the pole mass. In this 
way a precise evaluation of the total charmless semileptonic BF can be made directly using 
the on-shell scheme, where the u = 1/2 renormalon ambiguities of the pole mass and the 
series in Gu in Eq. (4.2) are regularized using the PV prescription. According to Sec. 2 in 
Ref. [52], this results in the following values^ ^: 

mfelp^ = 4.89 ± 0.05 GeV G}/^ 



0.928 ± 0.002, (4.3) 



corresponding to Nj = 4. Here the error in the pole mass is dominated by the error 
on the short distance mass m^^ = 4.20 it 0.04 (as quoted in Table 1). The uncertainty in 
Ttotai (-S — > Xuli>) /|Vub|^ is largely determined by this error. 

^^It is crucial to defined both mt and GI/^ using the same prescription; the product nib GI/^ is prescription 
independent. Throughout this paper we use the PV prescription when refering to these parameters, although 
this will not be explicitly reflected in the notation. 
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Gf = 1.16637 • IQ--' ctcm = 1/(130.3 ± 2.3) 

iKsl = (41.13 ± 0.63) • 10-3; l^tbl = 0.999119 ± 0.000026 [67] 

Mw = 80.388 ± 0.035 GeV 

mf^{mf^) = 165 ± 5 GeV 

as{Mz = 91.19) = 0.1176 ± 0.0020 =^ rj = as{Mw)/as{ij) ~ 0.561 

~ 0.15 GeV 
m^(mf^) = 1.295 ± 0.015 GeV [68] 
mf^(mf^) = 4.20 ± 0.04 GeV [68,69] 

Table 1: Values of input parameters in the calculation of the BF and their uncertainties. 

A useful tool in the calculation of the total B — > width is to normalize it by the 
semileptonic width; both b ^ c and b ^ u have been used in the past. In this way one 
exploits the better control one has on higher-order corrections in the semileptonic case, 
and reduces the sensitivity to the pole mass. There is however an important difference 
between the radiative decay width (4.1) and the semileptonic one (4.2) in the power of the 
pole mass. Thus, while the ratio 

V{B — > X,7, E > Eo)/ r,„,,, {B ^ XJU) (4.4) 

itself is an observable, and therefore renormalon free, it still contains an explicit dependence 
on the pole mass through {m^ / rn^°^^)'^ . This has been dealt with in the past by resorting 
to some alternative renormalon-frce mass scheme, a step that unavoidably introduces some 
uncertainty that is hard to quantify. Here we suggest an alternative procedure that utilizes 
the charmless semileptonic width to eliminate the explicit power dependence on the pole 
mass and yet does not require any additional mass scheme. This procedure is explained 
below. 



Using the semileptonic width 

In order to use the semileptonic result (4.2) to normalize the radiative one (4.1) in a renor- 
malon free manner, we first write the series for F over G^^ . By virtue of the cancallation 
of renormalon ambiguities in Eq. (4.1) and in Eq. (4.2), this ratio is renormalon free, so 

it can be evaluated at fixed order. To recover the result for the radiative decay width in 

3/5 

Eq. (4.1) we then multiply the series for F/Gu by the PV result [52] for G„ quoted in 
Eq. (4.3), namely gH^ = (0.928)3/5 and use the corresponding PV pole mass in the overall 
m\ power (the product m^G^^ is prescription-independent, see Eq. (4.2)): 



T{B Xsi,E > E^^) = \V,^VCf (mf (m,))' x 

\ ^ , (4.5) 

F{E^^)/G^J'> 



mlGT 



. Fixed Order 
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with 

^ /o(/^) + (fiit^) - ^^i/o(/^)) ^ (4.6) 

Fixed Order \ ^ / ^ 

+ (/2(m) - - ^^2(m)/o(a^) + §s? /o(m)) (^)' + • • • • 

Similarly, in order to use the resummed calculation of the normalized integrated spectra 
for individual matrix elements, namely Gjj(£'o)/Gy(-E'min)5 we define 

Fij {E^in) = Cf{^l)Cf{^i) Gij {Er^in, m) (4.7) 

= /o'(/^) + /r (Z^)^ + /2^(m) (^) ' + ■■■ 

TT \ TT / 

and write 

(4.8) 



Resummed 

This formula will be used for the numerical analysis that follows. In Fig. 5 above we show 
the contribution from each of the terms Fjj. Eq. (4.8) will allow us to study the theoretical 
uncertainties associated with different ingredients as a function of the cut. 

The total BF and the theoretical uncertainty 

Using Eq. (4.5) we obtain the following result for the total BF: 

BF{E^ > mb/20) = (357 ± 40(^) ± 19(^^) ± 21(p,,,„,)) • IQ-^, (4.9) 

where, as explained above, we used (4.6) to NNLO where /i(/x) is the full NLO coefficient 
while /2(/ti) is the /3o part of the NNLO corrections to the matrix elements, with the 
parameters of Table 1, which yield (4.3). The central value is based on = nib and //c = 
mh/2, and the three errors represent the variation of (1) /j, and (2) Hc as explained below, 
and (3) the parametric uncertainty in m^^ and ag^^ according to Table 1, respectively. 

Let us discuss now the theoretical uncertainty in some more detail. Having separated 
the calculation into that of Wilson coefficients on the one hand, and matrix elements on the 
other, the stability of the result with respect to the renormalization point of the operators 
becomes an essential measure of the accuracy of the calculation. It should be emphasized 
that by modifying /x in (4.1) or in (4.8) both the renormalization scale of the coupling and 
the factorization scale (i.e. the separation between the coefficient functions and the matrix 
elements) is changed. This necessarily implies reshuffling of contributions between different 
matrix elements according to the anomalous dimension matrix. 

It should be noted that there are various possibilities for choosing the renormalization 
point of the short distance mass in (4.8), which are of course reflected in the anomalous 
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Figure 6: The scale dependence of the total BF (the BF with > ■mb/20) at LO, NLO and 
NNLO, according to Eqs. (4.5) and (4.6). At NNLO only Poag^ terms are included in /2(m)- 



dimension matrix. We have chosen to fix this scale at the pole mass mb, rather than to 
vary it with fi. 

Figure 6 presents the variation of the total BF as a function of fi, with the default 
parameters of Table 1 and with ^/z = nicifJ^c) / T^b = 0.2215 (see below). We see that at 
leading order there is a large scale dependence (that of {Cj^{fi))'^) which get significantly 
smaller at NLO and NNLO. Note that at NNLO only /^qOs^ corrections are included here. 
The reason for suppressing the known O(as^) terms for G77 beyond this approximation, is 
that there is a large cancellation, at each order in the expansion, between G77 and G27- It 
is therefore essential to include the same type of corrections for both^'^. Note that there is 
no significant reduction of the scale dependence by including the PoOg'^ NNLO corrections. 
This is a reflection of the fact that important NNLO corrections are still missing: the 
leading scale dependence at each order, 0{as'^ln^ rrih/ fi), is associated with the evolution 
of the operators and their mixing. 

Another important source of uncertainty in the calculation of the total BF [7] is the 
renormalization point fj,c of the charm mass rric entering through the ratio z = m'^{fic)/^l 
into the expressions for the matrix element involving O2 (and Oi), notably G27. We 
have chosen the central value for z as ^/z = 0.2215, corresponding to an intermediate 
renormalization point, /ic — ?7ib/2 GeV, which is in between the b mass and the c mass. 
This is close to the central value chosen in [7]. For the error estimate we vary z in the 
range suggested in that paper, namely 0.18 < ^/z < 0.26. This uncertainty will reduce 
once the calculation of the matrix elements involving O2 is complete. 

Let us also note that the results of Ref. [49] for the /?o«s^ corrections involving charm 
loops were computed as an expansion m z = mc{^c) /^b- When using this expansion we 

^^For the default values of the parameters and with /i — mt, using the full G77 NNLO correction instead 
of its Poets'^ component (while keeping other matrix elements with the Pools"^ according to [49]) the BF 
increases by ~ 2.5%. 
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assume that it is valid at the physical value of z, an issue that was checked in Ref. [49]. 
Given that z is already varied over a large range, we do not consider the residual error 
associate with using a finite-order expansion in z as an independent source of uncertainty^^ . 

We note that the central value of Eq. (4.9) is somewhat lower than that of previous 
estimates. This is despite the fact that NNLO corrections increase the BF: using the same 
procedure at NLO we get BF(£:^ > mfe/20) = 322 • 10"^. These results can be compared 
with the result by Gambino and Misiak [7]: BF(£'-y > mf,/20) = 373 • 10~^, who used the 
NLO calculation. The main change is due to the different procedure for normalizing the BF 
by the semileptonic width. Note, in particular, that Ref. [7] evaluated the factor m^^^/m^"^*^ 
appearing in (4.4) by replacing it with mf^ /m\^ . At NLO this replacement has no effect 
on the expansion since the expansion oi / mJ^°^^ starts at ©(ct^^), however, numerically, 
this ratio is ~ 0.95. This explains much of the difference in the central values between 
the two estimates. It should also be noted that the uncertainty we find from varying the 
renormalization scale /n, ±11%, is larger than those in Refs. [7,8], and consequently the 
combined uncertainty on the total BF is about ±13.7%, somewhat larger than previous 
estimates. 

Finally, we point out that certain small corrections have been neglected in our cal- 
culation. The largest amongst these is the electroweak corrections whose effect has been 
estimated [7] to be —3.8% of the BF. Another correction is the power suppressed A^/m^ 
that was estimated as +2.5% of the BF [7]. As soon as more complete NNLO results 
become available these corrections should be taken into account as well. 

4.2 BF for > Eq 

In Fig. 7 we present the BF as a function of the cut, BF(i?^ > Eq). The uncertainty band in 
this figure as well as in Figs. 3, 8 and 10 to 12, indicates the theoretical uncertainty obtained 
by varying separately the following parameters, and summing the respective uncertainties 
in quadrature: 

• The renormalization scale of the operators (and the coupling) in the calculation of the 
total width according to Eqs. (4.5) and (4.6) between m},/2 and 2mi, where fi = 

is the default. 

• The renormalization point /ic of the charm mass rndp^c) entering the matrix elements 
involving O2 and Oi. As mentioned above we use: 0.18 < mdnc) < 0.26 m^. 

• The renormalization scale in the matching coefficients of the resummed spectra 
(A. 43), between mf, and m{,/2, where /x = rrib is the default. 

• The parameters controlling the details of the quark distribution function (corre- 
sponding to the leading "shape function") Bs{u), which enter the Sudakov factor 
in Eq. (2.14), with the corresponding power corrections, namely (NA/m,b)^ for any 
A; > 3. Specifically, as discussed in Sec. 4.3 below, we vary C3/2 in Eq. (2.9) (see 

^^See Note Added concerning a new publication on the subject that appeared upon completion of this 
work. 
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Eo [GeV] E„ [GeV] 

Figure 7: Left: The central value and uncertainty on the branching fraction ratio BF(i<^^ > Eq) 
as a function of the cut on the photon energy Eq. The black central line is obtained with the 
standard values of the parameters, while the (red) dot-dashed line is obtained by the same choice 
of parameters as for the dot-dashed line in Figs. 10 and 11. Right: The breakdown on sources to 
the uncertainty in the branching fraction. 



Eqs. (3.27) to (3.29) in Ref. [52]) and /^^ in Eq. (4.12) below within a range 
that is determined based on the support properties of the spectra, according to 
CutA, CutB < 0.01, where the functions Cut A and CutB are defined in (4.13). 

• The short distance parameters, af^ and m^^, within their error ranges as indicated 
in Table 1. 

Let us point out that the effect of the charm mass on the running of the coupling is 
neglected throughout the calculation of the spectra. We use Nf = 4, and ignore the charm 
threshold. The consequences of this approximation have not yet been studied. They 
certainly worthwhile investigating since at large the coupling in the Sudakov factor is 
effectively evaluated at scales below the charm mass. We further neglect non-perturbative 
effects that go beyond the summation of the leading powers, [N K/rri},)^ . This includes 
the so-called "subleading shape functions" , namely power corrections that are suppressed 
by 1/A compared to the ones we include, as well as the non-perturbative structure of 
the final state, which can be viewed as power corrections in (NA'^/mf), or directly in 
momentum space as exclusive states. 

Excluding very high cut values, the largest uncertainties shown in Fig. 7 have very 
little dependence on the cut. This includes, in particular, the renormalization scale of the 
operators and the value of mc{Hc)- The cut dependence of these contributions, is entirely 
due to the change their variation induces on the relative weight of the different matrix 
elements in Eq. (4.8). Furthermore, Fig. 7 shows that the uncertainty associated with the 
values of the short distance parameters, is also insensitive to the cut for sufficiently mild 
cuts. It is therefore useful to analyze separately the cut-dependent uncertainty. This is 
done in Fig. 8, where we display the partial BE ratio BF{E^ > Eq)/BF{E^ > mh/20) as 
a function of Eq. This ratio is needed, and can be used for extrapolating experimental 
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Eo [GeV] E„ [GeV] 

Figure 8: Left: The central value and uncertainty on the branching fraction ratio BF{Ej > 
Eq)/'BF{Ej > .2445GeV) and the branching fraction itself (insert) as a function of the cut on 
the photon energy Eq- Right: The breakdown on sources to the uncertainty in the branching 
fraction ratio. For low values of the cut, the uncertainty is dominated by the contribution from 
a variation in the renormalization scale in the matching coefficients (A. 43). This relatively large 
effect is originating from the matrix elements other than 6*77 that are known to NLO accuracy 
only. For Eq < 2.1 GeV, the parameters controlling the behavior in the Sudakov region contributes 
only little to the overall uncertainty. The (red) dot-dashed lines are obtained by the same choice 
of parameters as for the dot-dashed line in Figs. 10 and 11. 



Eo (GeV) 


default 


mill 


max 


1.6 


1.028 


1.025 


1.058 


1.7 


1.034 


1.031 


1.071 


1.8 


1.045 


1.041 


1.090 


1.9 


1.062 


1.056 


1.117 


2.0 


1.092 


1.083 


1.160 


2.1 


1.150 


1.134 


1.238 


2.2 


1.287 


1.247 


1.414 


2.3 


1.692 


1.556 


1.965 



Table 2: Extrapolation factors for the B — > X^i BF, m{E^ > mb/20)/BF{E^ > Eq) 



measurements from the region > Eq to the whole of phase space. Some numerical 
values for the extrapolation factor are presented in Table 2. As expected, the uncertainty 
of the extrapolation increases as the cut is increased. However, the rate of increase is not 
too high: precise experimental measurements with fairly stringent cuts, such as E^ > 2.0 
GeV, can still be useful. As shown in Fig. 8, up to cuts as high as Eq = 2.3 GeV, the largest 
source of uncertainty is still the dependence on the renormalization scale in the matching 
coefficients of Eq. (A. 43). At higher cuts the rrn, dependence and the details of the quark 
distribution function, Bs{u) and the associated power corrections become dominant. For 
the higher moments of the photon energy these ingredients become important already at 
milder cuts. We therefore discuss these issues in more detail in the next section. 
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'3/2 



Figure 9: The region in the u = 3/2 renormalon residue and power corrections plain, (C3/2, /^^), 
where the spectra conforms with the physical support properties. For example, in the most restricted 

region, marked by circles, both the requirement on the BF CutA(C3/2, ./^^) < 0.01 and the one on 
the average energy CutB(C3/2, /^^) < 0.01 are satisfied (CutA and CutB are defined in Eq. (4.13)). 



4.3 Spectral moments for > Eq 

Spectral moments defined over a restricted energy range E^^ > Eq [17], have proven to be a 
useful tool for comparison of the spectrum between experimental data and theory [12-14]. 
This comparison is important for a few reasons. First, it allows to test the extrapolation 
of B — > Xgj from the region of measurement. Second, more generally, it allows to test 
our understanding of the Sudakov region, which is a critical ingredient in extracting |Vub| 
from semileptonic decays, B — > Xylv. It that case the extrapolation from the region of 
measurement into the whole of phase space is very large, a factor of ~ 2 — 5 depending on 
the specific kinematic cuts applied, and thus the sensitivity to the details of the spectrum 
in the peak region is significant. 

Following Ref. [17] we consider here the average photon energy with a cut, namely 




(4.10) 
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Figure 10: Left: The prediction for Mb — {Ej)e-,>Eo as a function of Eq compared with data from 
Belle and BaBar. The green band indicates the uncertainty, and is obtained by adding contributions 
in quadrature. The default choice of parameters, corresponding to /i = ni^^ in Eq. (A. 43) and 
(C3/2,/^^) = (1,0) in Eqs. (2.9) and (4.12), is shown as a fuU (black) line. The dotted line is 
obtained by choosing fi — to^^/2. A slightly better description of the data is obtained upon 
choosing /i = m^^/2 and {C3/2, f^^) = (6.2, .3). This result is indicated by the (red) dot-dashed 
line. These choices all have equal theoretical justification. Right: The breakdown on sources for 
the uncertainty. 



and higher truncated moments, defined by: 



It is obvious that the higher the cut, or the moment n considered, the higher is the 
sensitivity to the fine details of the peak and, eventually, the endpoint region, ~ Mb/2. 
In contrast to the BF with mild cuts, where one can obtain a purely perturbative prediction, 
when considering the moments or high cuts one must take account of power corrections. In 
DGE, this can be judged by the sensitivity of the Borel integral to values of u away from 
the origin. In our parametrization of the soft function Bs{u) (see Eqs. (3.27) to (3.29) in 
Ref. [52]) this directly depends on Since Bs{u) at u ^ 3/2 is not known, and since 

the significance of power corrections in (2.10) directly depends on what it is assumed to 
be, it becomes obvious that these two aspects must be addressed together, and that an 
additional constraint would be needed. Here we propose to use the support properties, 
namely the vanishing of the spectrum for E.y > Mb/2 for this purpose. 

In taking power corrections in the peak region into account one has to find a com- 
promise that allows a sufficiently accurate description of the spectrum and yet involves a 
sufficiently small number of non-perturbative parameters. In general, starting with pertur- 
bation theory, the closer the region of interest to the endpoint, the harder it is to find such 
a compromise. Here we wish to explore the relation between C3/2 in the parametrization 
of Bs{u) entering Eq. (2.14), the power corrections and the support properties. As a first 
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Figure 11: Left: The prediction for the second central moment \{{E'y)E^>Eo ~ E-yY) ^ 

a function of Eq compared with data from Belle and BaBar. The band indicating the uncertainty 
and the various lines are obtained with the same choice of parameters as the corresponding ones in 
Fig. 10. Right: The breakdown on sources for the uncertainty. 



exploration of this issue, it is reasonable to consider a single non-perturbative parameter. 
On the other hand, since we would like to explore extreme cuts — the support properties 
— there is no justification in using only the leading power term, A; = 3 in Eq. (2.10). For 
~ Mb/2 such hierarchy does not exist. In order to use Eq. (2.10) we therefore take a 
sum of all renormalon ambiguities, all weighted by one non-perturbative parameter, /^^: 



Sud^ \N,mf,] 



Sud^ \N,mf,] 



X 

PV 

T{k/2) 



Bs{k/2) i?W(Ar,A:/2)|, 



(4.12) 



where /^^ is expected to be of order 1. Eq. (4.12) includes a sum of all powers of (NA/mi,)^ 
with k > 3, which a priori (depending on Bs{k/2)) may all be relevant for ~ Mb/2. 
Going over to milder cuts, where only the leading power is relevant, /^^ can be identified 
with in (2.10). 

We already know from previous studies that very large values of C3/2 or of power 
corrections cannot be accommodated with the support properties, E^ < Mb/2. Here 
we would like to translate this into a concrete constraint. To this end we require that 
neither the spectrum nor the first moment extend too far beyond the physical endpoint. 
Considering all possible values of C3/2 and /^^, we implement this requirement using the 
following cuts: 



CutA(C3/2,/P'') 
CutB(C3/2,/P^) 



T{B^Xs-/,E^ > Mb/2) 



T{B^Xs-f,E^ > Mb/20) 

^ {Ej)e^>Mb/2 



Mb/2 



(4.13) 
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Figure 12: Left: The prediction for the third central moment \{{E'y)E^>Eo ^£'7)'^) a 

function of Eq- The band indicating the uncertainty and the various lines are obtained with the 
same choice of parameters as the corresponding ones in Fig. 10. Right: The breakdown on sources 
for the uncertainty. This plot is terminated when the prediction for the third moment becomes 
consistent with zero. 



In Fig. 9 we have indicated by plus's the region in the plane of (C3/2, /^^) for which 
the spectra conform to CutA(C3/2, /^^) < 0.01 with the remaining parameters at their 
default value. We have also indicated by circles (asterisks) the region for which the spectra 
also obey CutB(C3/2, Z^^) < 0.01 (0.02). 

We observe that there is a reasonably large range in the (C3/2,/^^) parameter space 
that conforms with the physical support properties. The details depend on how stringent 
the constrains on (4.13) are. What is general, however, is that acceptable spectra have 
power corrections that are typically of the order of the renormalon ambiguity: as shown 
in Fig. 9 /^^ ^ 1, except at very small C3/2 where the power correction essentially has no 
effect. 

Having excluded too large contributions to the Sudakov factor from large u values 
and too large power corrections based on the support criterion, we still have a variety of 
spectra whose properties, i.e. the first few cut moments, are different. In the analysis of 
the theoretical uncertainty in this paper we allow (C3/2, /^^) to vary within the region of 
Cut A, CutB < 0.01. It is reassuring that the power corrections satisfying these require- 
ments are also of the order of the renormalon ambiguity, or smaller. Note, on the other 
hand, that for most of the range in C3/2 vanishing power corrections are excluded as well. 

In Figs. 10 to 12 we show the cut dependence of the first three central moments, 
Mb — {E^) ^ where the average energy is defined in (4.10), the variance n = 2 in 
(4.11) and the third moment n = 3 in (4.11), respectively. The bands in the figures on 
the left hand side represent the estimated theoretical uncertainty obtained by the same 
procedure used in Sec. 4.2. The various contributions to the uncertainty are presented in 
the plots on the right. As expected, the uncertainty associated with the behavior of the 
quark distribution function, Bs{u) at u ^ 3/2, and the corresponding power corrections 
increases as the cut is raised. It also gets more significant of course for higher moments. 
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In the average energy and variance figures, Figs. 10 and 11, we also present experi- 
mental data [12-14] by Belle and two analysis by BaBar. Fits to data, with simultaneous 
variation of mi,, C3/2 and f^^, can be very useful for the measurement of the bottom mass 
and for testing the theoretical description of the peak region. We do not perform any such 
fits^^ here, and the comparison with data is merely qualitative. 

The default line (full line) in Figs. 10 and 11, corresponding to the C3/2 = 1 choice 
made in Rcfs. [17, 52] with no power corrections /^^ = 0, seems to agree very well with 
the Belle data, and not as good with the BaBar data that often has smaller errors. In 
particular, in the average energy plot, the theoretical curve lies quite significantly above 
the BaBar data point at Eq = 2.26 GeV; in the variance plot it is about two sigma below 
the BaBar data points at Eq < 2.1 GeV. When making this comparison one should take 
great care in interpreting theoretical uncertainties, which do not reflect probability. For 
example, there is absolutely no theoretical preference to the choice C3/2 = 1 with /^^ = 
as compared, for example to C3/2 = 6.2 with /^^ = 0.3. Similarly, there is no preference 
for choosing the renormalization scale in (A.43) as ^ = nih (our default) as compared to 
H = mb/2. Short of further calculations, these choices remain arbitrary to a large extent. 
We note that with these different choices there is good agreement with BaBar data. This 
is demonstrated by the dot-dashed lines in Figs. 10 and 11. 

5. Conclusions 

We presented here a calculation of the B — > Xg^' branching fraction, as well as the first 
few spectral moments as a function of a cut on the photon energy, using the framework 
of Dressed Gluon Exponentiation. Building on our previous work [17,52], and on recent 
progress in fixed-order calculations [39,42,44,45,49], we now have accurate predictions 
for these observables and good understanding of the theoretical uncertainties and their 
dependence on the cut. 

We made progress on several different aspects of the calculation: 

• In matching the resummed spectrum into the fixed-order expansion (Appendix A), 
we made full use of the available NNLO results of the G77 spectrum [39,42]. In 
this sector wc therefore have a complete NNLO with NNLL accuracy. As shown 
in Fig. 1 the DGE spectrum does not vary much in going from NLO to NNLO, 
indicating that all important higher-order corrections are indeed resummed. The 
largest relative variation is along the tail of the distribution where NNLO corrections 
in the matching coefficient are important (see also Fig. 3). 

• We developed a method to perform Sudakov resummation without violating the an- 
alytic structure of the perturbative series in moment space, and therefore, without 
generating artifacts away from the Sudakov region. Using the Sudakov factor (2.14) 
the resummed spectrum can therefore be used down to small E^ where it matches 
into the fixed-order result (See Fig. 2). 

^^Fits should of course take into account the separation between statistical and systematic experimental 
errors (which are summed here in quadrature) as well as correlations between the data points. 
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• Matrix elements of other operators in the effective Weak Hamihonian and their in- 
terference with O7 are known in full to 0{as) only. However, we also know that 
independently of the nature of the short-distance interaction, all important contri- 
butions in the peak region, i.e. ones that do not vanish as 0{1/N), necessarily involve 
the same Sudakov factor. Here we computed the spectra of individual matrix ele- 
ments Gij assuming the same Sudakov factor multiplied by 0{as) hard coefficient 
functions that depend on the operators. The resulting effect on the spectrum is shown 
in Fig. 5. 

• Wc performed a first numerical study of the relation between renormalon contribu- 
tions to the quark distribution function, power corrections, and support properties. 
To this end we used the parametrization of the soft anomalous dimension function 
Bs{u) proposed in Ref. [52] (see Eqs. (3.27) to (3.29) there), which is consistent with 
all available constraints for u ^ 0, and at u = 1/2 and u = 1; it includes a sin- 
gle parameter C3/2 that controls the u ^ 3/2 region, which is not well constrained 
theoretically. Wc also used a simple formula with a single parameter for the cor- 
responding power corrections, Eq. (4.12), which nevertheless takes into account all 
powers of {NA/nib)'' with k > 3. Under these assumptions we showed that there 
exists a range in the parameter space where the computed spectrum conforms with 
the physical support properties. We observed that the size of the power corrections 
is of the order of the renormalon ambiguity, as expected on general grounds. 

• We proposed a new method to evaluate the total B — > Xg^ width that utilizes 
the semilcptonic width, where the leading renormalon ambiguity has been explicitly 
dealt with using Borel summation [52]. Our new method avoids using any addi- 
tional mass scheme. The result for the total BF is summarized in Eq. (4.9). It is 
consistent with previous determinations [7, 8] . Unfortunately, despite having partial 
NNLO information on the matrix elements [44,45,49], the theoretical uncertainty is 
still large, ±13.7%. We observed large cancellations between G77 and G27 contri- 
butions, which vary depending on the renormalization scale. Therefore, significant 
improvement is expected upon completion of the NNLO calculation of G27. 

• Finally, we devised a method to compute the BF with a cut on the photon energy, by 
combining the resummed calculation of the spectrum of individual matrix elements 
with the proper weight, given by the fixed-order calculation of the total BF, Eq. (4.8). 
This framework facilitates the analysis of theoretical uncertainties associated with 
different ingredients, which are known at different orders, and have different role 
depending on how stringent the cut is. 

It is of theoretical as well as practical interest to understand the behavior of higher- 
order corrections. The overwhelming dominance of running-coupling corrections at ©(a^^), 
and the very late settling of the leading logarithmic behavior at large x [39] , might be inter- 
preted as a signal that Sudakov resummation is irrelevant. In order to address the relative 
significance of higher-order corrections of different origin wc presented in Sec. 2.3 the all- 
order resummation of the G77 spectrum in the large-/3o limit, based on the calculation of 
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Ref. [18]. When working directly in momentum space (x space) there are no renormalon 
ambiguities in the real-emission result, but there are convergence constraints on the Borel 
integral. Consequently, the Borel sum does not exist for GeV. As shown in Fig. 1 

and 4, the large-/3o sum does not give a viable description of the spectrum in the peak 
region. This stands in sharp contrast with DGE, which through real-virtual cancellation 
and exponentiation in moments space, accounts for multiple soft and collinear radiation. 
Thus, despite the late settling of the leading logarithmic behavior at large x up to 0{as'^), 
such logarithms are important at higher orders. It is the combination of soft gluon resum- 
mation with the resummation of running-coupling corrections that open the way for the 
quantitative description of inclusive decay spectra. 

Finally, the results of the present paper can help in various ways in understanding 
and using the B factory B — > Xg'y data. In order to make comparison between the BF 
measurements with a given kinematic cut, > Eq, and the theoretical result of Eq. (4.9) 
we provided in Table 2 estimates of the extrapolation factor as a function of Eq. This 
information, as well as all other details of the spectrum and the moments computed here, 
can be obtained using a C++ program that is made publicly available^^. Beyond the issue 
of the BF itself, combining the theoretical calculations presented here with B — > Xg^ data 
is extremely valuable for other aspects of flavor physics: it allows a precise determination 
of the the b-quark mass, as well as testing and improving the description of inclusive decay 
spectra in the Sudakov region by quantifying Fermi-motion corrections. This is particularly 
important for the determination of |Vub| from the B factory measurements of charmless 
semileptonic decays [15,74]. 

Note Added: Upon completion of this paper there appeared two new papers addressing 
the calculation of the total BF [75,76]. In particular, Ref. [76] estimates the 0{al) contri- 
bution of the four-quark operator matrix elements using extrapolation from large rric- This 
goes beyond the large-/3o results of Ref. [49] that we used here. The final result for the 
total BF in [75, 76] is consistent with ours, but has a lower central value and a significantly 
smaller Tincertainty from varying the renormalization scale. The reasons for that are not 
yet clear. 



A. NNLO matching of the resummed G77 spectrum 

A.l Expansion of the Sudakov exponent and a basic NNLO matching formula 

Let us begin by recalling the O(as^) result for the Sudakov exponent [17]. Expanding 
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Eq. (2.6) to this order one obtains Eq. (2.5) with 



Ei{N) = -y ( 5?(7V) - S2{N) ) - (61 - di)Si(iV) 



(A.l) 
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and where the coefficients^^ of the Sudakov anomalous dimensions, defined in the MS 
scheme, are [17,39,43,53-55]: 
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In the following we will compute the matching coefficient functions H{as{mi,), N) and 
AR{as{mi)),x) in Eq. (2.4) to ©(a^^) based on the known NNLO expansion [39] and the 
expansion of the exponent in Eq. (2.5) with (A.l). 

To this end we will need the color decomposition of the NNLO result in Eq. (2.1) 
(note that an overall factor of Ci;' was extracted there). The NLO coefficients entering 
R{as{mb), x) and V{as{nni)) are ri{x), which is is given explicitly Eq. (B.4) below, and 
ki = —31/12. The NNLO ones will be decomposed as follows: 



k2 = Nfk'^f + Cpk^^ + CAk^^; 



N. 



rf '(a;) = Nfr^' [x) + CFr^"'^'°^(x) + Cav^ 



C Fusing. 



CA,sing. 



(^); 



(A.4) 



r'^'^ix) = Ar^rf^"'*(x) + Cpr^^'^^ix) + C^r^^ (x). 



^^'Note the diflFerent normalization compared to the standard harmonic sum. 

^^For the precise relation with Bs{u) and Bj{u) see Ref. [17]. Note that the notation used there in Eqs. 

n-l. 



(2.5) through (2.7) is: ^„ = C'fPS a„ etc. 
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The explicit expression for r^^'"'^ (x) = —^r2'^'''''^'{x) is given in Eq. (B.5) while the expres- 
sions for the other terms, corresponding to the Cp and CfCa color factors, can be found 
in Eqs. (7) and (8) in Ref. [39], respectively. Below we quote only their singular parts, 
r2°^ (a;), which were derived already in [17]. The regular parts, r2^ {x), are used here but 
owing to their length we do not write their explicit expressions here; these can read off Eqs. 
(7) and (8) in Ref. [39] by removing the 5 and the plus-distribution terms. Using this color 
decomposition one finds the constants kn by integrating r'j^^ {x) according to Eq. (2.2) [39]: 

k^f = ^ - k + — = 1.237; k^^ = -4.795; A;^^ = 1.216. (A.5) 
^ 48 3 16 '2 '2 V ; 

Having set the notation let us consider now the partially-integrated G77 matrix element 
at O(as^), as obtained from Eq. (2.3) with (2.1): 
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Defining the moments of r^'°^ (a;) according to Eq. (1.2), 
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Thus, we can express Eq. (A.6) as an inverse-Mellin transform: 
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Next, wc can incorporate the resummation of Sudakov logarithms by rewriting Eq. (A. 9) 
in terms of the Sudakov factor of Eq. (2.6) according to the general form in Eq. (2.4): 



G77{Eo,mb) 
_ 077(0, nib) 



Resummed 

X exp 



1 

2TTi 



C+JOO 



dN 
N - 1 



2Eo 
nib 



l-N 



Sud{N,mb) 



(A.IO) 



CFk^^+CF{h-CFkl/2)^^^ +. 



+Cf 



f 

Jxn 



xo=2Eo/mb 



dx r"*' (x) 



Ois{nib) 



TT 



+ ri'ix) 



a 



TT / 



+ 



- 41 - 



Note that f2°'{x) introduced here is different from r2°^ (a;) used in Eq. (A. 9) since the 
exponentiation generates some additional integrable terms (see below). Such terms need 
to be subtracted-*^^ of r2^' (x) to avoid double counting. Of course, k2 is defined accordingly, 
kn = - dxf'^^ix). 

Upon expansion, this formula together with the coefficients of the Sudakov exponent 
was used to predict [17] the terms that diverge at N — > oo at 0{a1), see Appendix A 
there. Explicit calculations later confirmed these results [39,42]. Expanding Eq. (A. 10) 
using Eq. (2.5) one gets, under the inverse-Mellin integral: 
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Ef{N) + kiEi{N) Cf + E2{N). (A. 12) 



Eq. (A.ll) should be compared with the contents of the curly brackets in Eq. (A. 9). The 
0{as) terms are clearly identical, but the 0{a^g) terms are not. They are equal as far as 
the terms that diverge at — > 00 are concerned, yet they differ by 0{1/N) corrections 
proportional to C|, that are generated by exponentiating Ei{N). To account for this 
difference we have introduced the bar notation. Specifically ^2"^ (A^) is different from 
the originally defined i22°^ (A^) of Eq. (A.7) by some additional terms that are finite at 
N — s- 00: 
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The constant at N — > 00 is retained at its computed value by requiring ^2°^ (A) + ^2 = 
^sing.^^-j + A;2 + C(l/A), see Eq. (A. 20) below. In the remainder of this section we derive 
explicit expressions for the barred matching coefficients in Eq. (A. 10). 

Starting with Eq. (A. 12) and using Eq. (A.l) we get the following exphcit expressions: 



RriN) = - 2 ( St{N) - S2{N) ) - (61 - di)Si{N) 



lst{N)+^-{h-d^)Sf{N) + 



lk,-\s2{N)+^-{d,-hr]sf{N) 



+ (di - 61) ( ^S2iN) + ki) S,{N) + ^Sl{N) + ^kiS2 



1 



1 



- 7;Sm) +(di--h- -a2 St{N) + -S2{N) + d2-b2] -Si (AT) 



+ (-di + ^61 + ^02] S2iN) - ^SsiN) 



(A.14) 



^This is done in Eq. (A. 19) below. 



-42 - 



where we substituted ai = 1; other coefficients can be read off Eq. (A. 3). The additional 
terms contained in ^2"^ (-^) but not in i?2"^'(-/V), (finite terms at — > 00) are: 







dx (x^-i - 1) 



sing., integrable 



(x) = Cf 



+ (61 - di) (^(^S2iN) + Si{N) - ^SsiN) 
- i {S2iN) - Sf{N)) - Sl{N)S2iN) + {Ss{N) - 2(3) S,{N) ^ ^S,{N) + S^iN) 
+ Ulsl{N)-^S,{N) + ^S2{N) 



(A.15) 



In X space this corresponds to 

1 — a; 1 — X 

where the original terms, appearing under the plus prescription in Eq. (2.1), are^^ 

1 



(A.16) 



r-;'"^(a;) = Cf 



(:3-h{di-h) + —{bi-di) 

b J 1 — X 



+ { (di - hf -^-h 



(62 



6 

1 



Inl-x 3,, ,,ln2l-.T lln-Vl-x 

— ^ + - (di - 61) ^ -\ ^ 

1-x 2^ ^ 1-x 2 1-x 



1-x 



, ,ln(l-x) 31n2(l-x) 
+ (2di -61 -02)— ^ ^+ ^ ^ 



l-x 2 1-x 



(A.17) 



while the additional, integrable terms that arise from CpEf{N) are: 



sing., integrable 



(x) = Cf 



TT^ _ 1 

12 2 



{di - hif 



1-x 



1-x 



lln(x)ln2(l-x) 1 1 



1-x 



+ - 



TT 



2 ^ _ ^ ln(x) ( Li2(x) + — ) + 2Li3(x) - 2(3 



(A.18) 



It is now straightforward to obtain the modified (barred) r2^' (x) that enters the match- 
ing formula of Eq. (A. 10): 



f'i^\x) = r'^^{x)-r- 



sing., integrable 



(x). 



(A.19) 



^These terms are the same as in Eq. (A. 4) in [17], as they must be. 
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Note that only the Cp term is modified. Finally, the corresponding change in k2 can be 
easily determined by computing 

k2-k2= C dxrf (x) = - lim « (iV) - i?f « (A^)1 

^° /r"' (A.20) 

c 

Using Eq. (A. 5) we therefore find k^" = 1.7662, so 

h = -4.795 Ca + 1.237 Nf + 1.7662 Cp. (A.21) 
A. 2 Including more in moment space 

One advantage of the DGE, which is tightly connected to the fact that the calculation 
is done in moment space, is that the resummed spectrum smoothly extends beyond the 
perturbative endpoint x = 1 and tends to zero at a; = 1 + 0{A/mh), even in the absence 
of power corrections. 

For the differential spectrum to be smooth at x = 1, not just integrable, one must take 
into account in moment space not only terms that diverge (or are finite) for N — > oo, but 
also 0{\/N) corrections, corresponding in particular to powers of ln(l — a;). In the basic 
matching formula of Eq. (A. 10) these correction are still part of r^^^^x). In the following 
we rearrange the split between moment space and x space to incorporate all the terms 
that are finite for x — ^ 1 in moment space, leaving only subleading terms, that vanish 
at X — ^ 1 in a; space. In addition we include in moments space all running-coupling, 
O{l3oal) effects. 

One way to ensure this is, of course, full moment-space matching. One defines: 

R'^^ {N)= [ cix(x^-^ - l)f7 (x) = / dxx^-^r^^ {x) + kn. (A.22) 
Jo Jo 

and then 

GMEo,m,)^^r-^-m r^y-%,d(iV,n.,)x exp|c.i.-(iV)^ 



+ Cf 



1 M / _ N \ 2 

R^'^iN) -Cf{- {K{' iN)f + (R^^ iN) - ki)Ei{N) 



as{mb)\ ^ 



(A.23) 



Here we choose to implement full moment-space matching at NLO while splitting the 
NNLO corrections as follows: 

where f'^^-''"^^^'"^'^"^^(^x) vanishes as (1 — x) at x — > 1 (up to logarithms) and 

j^rcg., leading ^ f (^3; f-S- - I-^ing ^^^j ^^7V- 1 _ _ /" dx f'^^' ' '"'"''""'"'^x) 

Jo Jo 

= kn + f dxf''^''"""''''^{x)x^-^ (A.25) 
^0 
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such that at iV — > oo: K"" ''°"'''""(Af) — > while at finite N Rn""- (N) contains ah 
the terms that faU as (including 1/N times a power of In A'^). We shall not make this 
split for r"^' nor for the /?o contribution to f2^' (the so-called BLM term) — these will be 
fully contained in the moment-space expression. 

For the leading terms in fg"^' at x — > 1, Eq. (A.24) we find: 



,A.„.e,., leading ^^-j ^ _ 7 _ ^ ^ 13 ^^^^ _ ^) _ ^ 1^2(1 _ 



24 36 36 



(A.26) 



-C..e..,eaa..(^) = 7^"^ " " ^ + (- ^ - ln(l - x) + ln2(l - X) + i ln3(l 

eg., leading,^, _ ^ ^ _ f_ ^ ( ^ f_] l^^J - x) + - lli^l - x) . 



X) 



4 36 



9 



In moment space, the coefficients of the different color factors of FQ^ ' "'^ ""^{N) read: 



„iVy , reg., leading 



R. 



Cjr 5 reg. , leading 



R'. 



tjCa , reg. , leading 
K2 



(AT) = fef^ + ( - J- - ^ J ;,o(iV) + ^//l(iV) - ^M2(iV) 



7 TT^ 



13 



(A.27) 



12 4 



+ Ai2(iV) + -^3(iV) 



(iV) = ^2''- + (^C3 + ^-|^)Mo(iV) + 



-f + y) Mi(iV) + ^M2(iV), 



//,(iV) = /Jdxx^-i W(l-x), e.g. ^lo{N) = 1/iV, Mi(iV) = -{^,{N)+^e)/N etc. 

Making this split, and shifting the renormalization scale of the coupling in the matching 
coefficient to an arbitrary scale /x, we obtain: 



GT7{Eo,mb) 

_ G77(0,m;,) 



Resummed 



TT^ / T-r — 7 — - Sud(N, nib) X 

2 



TT 



/?o In ( ) ii;;^^ (Af) + i^-^^-'— (AT) 



-Cf[\ {R--{N)f + (R^iN) - h) E^{N) 



TT 



+ 



+ Cf I dr ^-S '-i'i-din, 

J xo=2Eo/mf, 



dx f-S.,sn.leadmg^^^ 



as{ix) 



TT 



+ 



(A.28) 



where the explicit expression for ^^^ {N) is given in Eq. (B.14), /3o = — \Nf and 

(A.29) 



_^reg. leading ^^^j = (i^R^ ^' '""'' {N) 



BLM 



+ CFii:^^''''«'''^"'''"«(Ar) + Ca[R^ 



Ca, reg., leading^^^ _|_ 11 ^A^/ , reg. , leading j 



non-BLM 
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where i?7'^^^ (iV) is given explicitly in Eq. (B.16), while the leading terms of the non-BLM 
color factors are given in Eq. (A. 27). As announced above, terms that are excluded from 
Eq. (A. 29) are 0{1/N'^) and do not involve running-coupling effects. These residual terms 
are included in x space through 

-reg., sublcading ^^-^ _ ^^-C'i? , reg. , sublcading ^^-j 

+ Ca (^ff l"«-l'"S(2;) + ^^^/. heading 

As usual, the renormalization-scalc dependence in Eq. (A. 28) can serve as a measure of 
subleading perturbative corrections at C'(af) and beyond. In the numerical analysis we 
will vary fi between mi,/ 2 and mb as one of several means to estimate the theoretical 
uncertainty. 

A. 3 NNLO matching under constraints on the analytic structure (J ^ 0) 

Let us describe now the matching procedure with the Sudakov factor of Eq. (2.14). The 

matching cocfhcicnts, just like the exponent, are constructed under a constraint on the 
analytic structure in moment space: no poles should appear for N > —J, and so the small-a; 
asymptotic behavior would coincide with that of the fixed-order result, dT/dx ~ x"^ . 
The expansion of the J-modified Sudakov exponent, Eq. (2.14), takes the form: 



Sud^ \N,mh) = exp < Cp 



TT \ vr / 



(A.31) 



with 

En{N) = En{N + J)- En{J + 1) (A.32) 

where for n = 1, 2 are given explicitly in Eq. (A.l). 

We now proceed to match Eq. (2.14) into the NNLO result. First, let us write a basic 
matching formula, in analogy with Eq. (A. 10), 

G77(-Bo,m6) 1 dN f2EoY-^ 



G77{0,mb) 



= _L r- (?^) sVd'^V.m,) (A.33) 



exp \c,h^ + C, {k, - Cklli) (^) 



X exp < CVfci + Cf I k2 - CfKI2 ) | ^ + • • 

dx (fT (x) + vT {x) ( ' + • ■ 

Jxo=2Eo/mb \ ^ \ 1" / 

Defining i?^"^' (N) as the coefficients arising from the expansion of the Sudakov factor times 
the matching coefficient in Eq. (A.33): 

5p(Ar) = Ei{N), 5f «(Ar) = (^EfiN) + hEi{N)\ Cp + E2{N) (A.34) 
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we can compute the modified regular parts as follows: 



f;^-s- (x) = r'^^- (x) - Ar'^""^- (x) (A.35) 
where Ar^'"^- {x) = f^°^- (x) — r^"^- (x) and where 

^sing.^^-j = f {x^-^ - 1) f^'^ lx), (A.36) 
similarly to Eq. (A. 7). It is straightforward to compute ki: 
ki-ki= f da;Arp (x) = - lim fe""' (iV) - ii^' ( A^) 

Jo N- — >oo L J 

= - Q^i'll + J) + {hi - di)Si{l + J) - 1^2(1 + J)) Cf 

(A.37) 

Because of the structure of En in Eq. (A.32) one can express R'^^^N) as^° 

Rt^\N) = Fn{N + J,J)- Fn{l + J,J)= C dxx^'+J-^ [fnAx)\+ - t dxx^ [fn,j{^)]+ 

Jo Jo 

= C dx {x''-^-l)x'fn,j{,^). 

Jo 

(A.38) 

This means the inverse McUin transform of i2^J"^ (A^) is readily obtained by multiplying the 
inverse Mellin transform of Fn{N,J) by x'^ under the plus prescription. Moreover, since 
for J = Eq. (A.38) must coincide with Eq. (A. 12) we have: 

Fi (AT, 0) = Rf''^- {N) ■ F2 {N, 0) = (AT) , (A.39) 

and therefore Eq. (A.38) implies^^ 

^p(iV) = i?P(Ar+J)-i?p(l + J); ^p(Ar) = «(A^ + J)-i?f «(1 + J) (A.40) 

where the explicit expressions of the functions on the r.h.s. are given in Eq. (A. 14) and 

Arr^ ix) = (x-^- l)rp (x) 

^ ^ (A 41) 

Ar^f^ix) = x-^ ff ^ (x) - rt^ {x) = x^ ^^j _ ^=ing.(^) 



^"As we show below fn,j{x) are in fact J independent upon substituting for fei in terms of fci. They are 
given by h,j{x) = ri'"'^ (x) and /2,j(x) = f2"*^ (a;). 

^^One can expUcitly compute iif 2^ (iV) and therefore Fn{N + J, J) using Eq. (A. 34). It is straightforward 
to verify that upon substituting for ki in terms of fci using Eq. (A.37), the dependence of Fn{N + J, J) on 
J appears only through N + J and the result coincides with Eq. (A. 14). 
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where r^'"'^ (x) and r2"*^ (x) and ^='"s.'mtogra c^^^ given in Eqs (A. 16), (A. 17) and (A. 18), 
respectively. Finally we find k2 by 



-k2= f dx^rf^ {x) = - lim i2f ( AT) - ( AT) 
Jo ^ — 



Cf 



ls^\l + j)-^-{d,-h,)S,\l + J) 



+ [\{di- hi) 52(1 + J) + (di - 6i) 5i(l + J) + i 52^(1 + J) + \ki S2{1 + J) 



+ /?0 



i51(l+J)3+L-l- 1)51^(1+ J) 



+ Q 52(1 + J) + d2 - 62) 5i(l + J) + (^-di + y + y) ^2(1 + - ^ 53(1 + J) 

(A.42) ' 

where we used Eq. (A. 37) and wrote ki in terms of ki. ^2 on the l.h.s is given in Eq. (A. 21). 

Having determined all the ingredients in the matching formula Eq. (A. 33), we can 
easily convert it to the preferred form where all the NLO terms and the leading NNLO 
terms are evaluated in moment space, in analogy with Eq. (A. 28): 



G77{EQ,mb) 



077(0, rrib) 



1 dN f2Ec^^-^ 



Resummed 



2 



exp <^ CpRr iN)^^ + Cf 



Po In 



mi 



Cf (\ (i?r (^))' + (i?r (^) - ^1) EiiN) 



TT 



+ 



_j_ (J^ I ^rn ;r''''S- ' subleading 



dx f-S-.-"'— 



as (a*) 



TT 



+ 



(A.43) 



where we used Eq. (A.34). Here R"^ {N) is 

R{^ {N) = f dx{x^-^ - iy{^ {x) = R{^ {N) + i?p (A) - ^p (iV) 
Jo 



i?P (A + J) - i?P (l + J) - Rf^' iN) 



(A.44) 



where the second line in based on Eq. (A.41), and i?^'"^ (A/") is given in Eq. (A. 40). An 
explicit expression for the regular part in x space can be obtained using Eqs. (A.35) and 
(A.41): 



f-'{x) = r-^{x) 



J sing., integrable 



(x) + (x-^-l)rr(x) 



(A.45) 
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Finally, ' '-^ing^^^ defined in analogy with i^-s-.'-ding^^^ ^^ 24) and (A.25). 

To this end we decompose the regular part: 



such that f-g-.-bieadingj.^-j vanishes as (1 - x) for x — ^ 1 and 

''''%x) (x^-^ - 1) - 

Jo 



(A.46) 



ding^^-j 



^^~reg.,leadi, 



dxf. 



"res., subleadinj 



(x) 



k2 + dxf 2 
'o 



reg. , leading / \ — 1 



(A.47) 



In contrast with r2''^'""^''''{x), however, we require that ■f2'''^'""^'''^{x) — and therefore 
also f-^^'-bieading^^^ _ ^.^ behave as x-' at x — > 0. To compute f2''«''''"'''°«(x) we 
therefore first extract x'^ out of 7^2'^' (x) of Eq. (A. 45) before expanding near x — > 1. Using 
x-^ = 1- J{l-x) + 0((1 - xf) and Eq. (A.19) we obtain: 



Finally, returning to Eq. (A.47) and using Eq. (A.25) we obtain: 

"'"^(AT) = k2 + j\x (^'^^■'"'^'''\x) + J{l-x)rf^ {x)^x^^^-^ 



^reg., leading 



(A.48) 



(A.49) 



N+J-l 



= k2-k2 + i2-^ ''-^-^(iV + J) + J / dx (l-x) r^^ix) x 

Jo 

where the explicit expression for ^^"^s-iieading^^^ given in Eq. (A. 29) with Eq. (A. 27). 

B. The normalized G77 spectrum in the large-/3o Umit: results 
B.l Expansion coefficients in x space 

The coefficients r^" of Eq. (2.19) are obtained upon expanding the Borel function in 
Eq. (2.17): 



-. 00 

-B{x,u) = J]r*(x)u"/n!; r*(x) = r*'^'"«(x) + r*'-*(x), (B.l) 



n=0 



To this end we need the expansion of the hypergeometric function, which is (see [66]): 



2Fi([l,l],[2-ii],x) = i^<j -ln(l-x) + n 



1 



ln^(l -x) + Li2(x) 



1 



(B.2) 



-S'i,2(x) - ln(l - x)Li2(x) + Li3(x) - ^ ln3(l - x) 
- S2,2{x) - ln(l - x)Li3(x) + ln(l - x)S'i,2(x) 



+^ ln2(l - x)Li2(x) + ln^(l - x) + 5i,3(x) + Li4(x) 



+ 
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where Nielsen integrals are defined by 



Sa,b{x) = 



ia-l)\b\ Jo e 



^ lln-^-^lOln^ll-Cx). (B.3) 



The resulting coefficients are as follows: 



1 1-x 4{l-x) 

rrix) = ln(l - ^) - \ + \ + \ 



/3o,.in,.,._^ 3hi2(l-x) , 13 1n(l-x) , (t:^ 85 



""2 ' ^^^"2 1-x +12 1-x 2471 

„/3oreg./ X _/2; l 1 \t-/\,/3x 3\ 2 



(B.4) 



, 25x 3 1 , , 19x2 49 55x tt^ 
+ — — - - 77: + TT ln(l - x) — + — + 



12 2x 12 2 y ' ' 12 24 24 6 (1-x) 

(B.5) 



(x) 
(x) 



7 ln^(l -x) 1 ln2(l -x) 
"3 1-x + 4 1-x + 



X + 1 + 



Li3(x) 



275 _ 
"36" " 

- 2x 



TT^X ln(l — x) 



+ 



29 TT^ 

36 



581 ' 
'72 



X 



1 X 

ln(l - x) + + 



3 10 
X 6 



3 ■ 3(l-x) 

Li,(x) + (-X - 1 - ^) 5i,2(x) + (-1 - ^) ln3(l - X) 



+ 
+ 

+ 



23x 9 13 x^ \ , 2/ 



19 x^ 7 TT^ 247 

\ 1 H 

6 X 6 36 

7r2 203\ 2 27r2 

"e" ~ ^ ~ 9(1 -x) 



+ 



407 TT^ 
"36" "^ Y 
923 

^ ~ 12 



27r2 



X 



3 (1-x 

139 77r2 
~~72 12" 



ln(l — x) 



(B.6) 
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A=ing.(^) _ 15 ln^(l - x) _ 35 ln3(l - x) ^(^_ 105 ^^ - x) 



1 — X 



12 l-x 



+ 



J 1 — X 



/6029 297r2\ ln(l - x) / tt^ 235 2 9557 
+ (vW"^J l-x +r30 + 72"'' "l32" + ^^^y l-x 



Jo, 



(x) = 3 1 + X + 



1 — X 



Li4(x) + 



1 — X 



Li3(x) + 



6 1 + X + 



1 — X 



2 \ 9 
6 ( 1 + X + ln(l - x) + - + 11 X - 10 

1 — X / X 



ln2(l - x) 



197 TT^ /121 7r2\ 21]^ , , 

-^-y + (^-yJ^+xJ^^^(^) 

+ 6 ( 1 + X + ) ln(l - x) — + 10 - 11 X - - 

1 — X/ 1 — X X 



+ 3 1+ x + 



l-x 



(5i,3(a:)-52,2(x)) 



15 15 X 



+ ^ + 



119 157x 21 



12 



^) ln3(l-x) 



/ 37r2 27r2 63 19x2 345 /297 Svr^X \ , o,, 

+ \ h X In^fl -x)+ 

V 4 l-x 2x 4 8V8 4;;^ ' 



12163 37r2- 57 47r2 



ln(l - x) + 



2x 3(l-x) 
4955 19 TT^ ^ 



203 tt" 



x^ + 



21 



12 2 ; V 12 

15715 / 55 TT^ 30227\ 



25 TT^ _ 14063 \ 7r_ 
216";' 12 



216 12 J 432 
1 -257r2 + 270C3-37r'' 497r2 
90 



+ 



24 



+ 



432 J 



l-x 



24 
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B.2 Expansion coefficients in moment space 



Defining the moments, as in (1.2), by 



M 



FT, Or 
N 



1 dror^N_i 



dx 



(B.8) 



and using Eq. (2.15) with the final expression in Eq. (2.17) we get: 



(B.9) 
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with 



B{N,u) 



(B.IO) 



5„ sinTT-u 



X {1 - u){2 - u) 



+ 



-(1 -Ax + x'^ 



1 1 

+ 



1 — x 1 — n 



+ 



2(1 - xf 

(l-«)2 



Fi [l,l],[2-n],x 



+(l-4a; + x^)- + 



1 (x + l)(a;2 - 3x + 1) (2-a;)x 2(1 - x) 



l-x (l-u)(l-a;) (2-«) (l-«)2 



The perturbative series to leading order in the flavor expansion takes the form: 



TT 



as{ml 



TT 



+ 



(B.ll) 



where at each order k we separate the coefficients into singular and regular parts, as in 
Eq. (B.l), namely 



RkiN) = Rt--{N) + R--iN). 



(B.12) 



Here both parts of the coefficients are defined with vanishing first moments at each order: 
i?p (7V = 1) = and R'^^ {N = 1) = 0. 
At NLO the coefficients are: 

i?p(iV) = l^i(iV)-^-i(M'(iV)+7^;)2 + ^ {nN) + lE) (B.13) 



LL 



NLL 



(B.14) 



9 

+^ + 



+ 



4N 2{N + iy 2{N + 2) 4 (A^ + 1) 2Ar2 



At NNLO we have: 



R. 



N-f ,sing. 



(iV) = 



-i(^(iV) + 7^)3 + U^{N) - ^\ (*(A) + je) - liMN) + 2C3) 



LL 



13 
+ 72 



+ {i^-^)(W+7^) (B.15) 



NLL 



NNLL 
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and 



oo 



mm) 



12 ^ r(iV + A; + 2) 

k=l 

1 13 31 1 1 / 1 



1 + 2N N(N + 1] 
3 + 3^ + 2- ^ ' 



fc2 



49 
+ 48 



+ 



1 



12A^3 72iv2 48 iV 8 \N N + 1 
11 7 1 11 



(B.16) 

19 1 11 
72 iV ~ 6iV2 



+ ^*i(iV)-^ ^ 



4Ar- 1 



(*(iV)+7E) 



12 AT + 2 72Ar + l 6 (A^ + 1)2 6 
24VA^ + 1 iVy ^ Vl6 144 AT 144 AT + 1 ^ 12 ^ S^'^ 



13 



1 



+ 



1 



+ 



1 



1 



+ 



7 



72 (AT + 1)2 12(A^ + 2)2 36A^ + 2 12(Ar + 1)3 144(A^ + 1) ' 
B.3 The small— X limit 

Upon expanding Eq. (2.3) in Ref. [18], or Eq. (2.17) above, in powers of x one finds: 
1 dTor 



dx 



large /3o 



^rduTiu)(^\ 1 5(^')(u)x3 + 5(^')(«)a;^ + 0(x5)) (B.17) 
2po Jo V»^6/ 



2tt 



35 agdn 

^ + i2ir + 



, , 5 ^ 1099a,/3o ^ 

" + '4+ Mo-^ + - 



x^ + 0(x5) 



with 



B^^ >(u) = es" 

TTU 3 



63' 



1 — u A — u 

5 1 4 
+ 



6 



7 



5 1 



21-'u 2-u 3-u 4-u 25-'u 



(B.18) 



Performing the Borel integration in Eq. (B.17) (with the default values = 4.875 GeV 
A = 0.332 GeV and Nf = 4) we obtain the following numerical coefficients of the small x 
expansion: 



1 dr 



07 



Vq-! dx 



large Pq 



0.0729 ar^ + 0.105 a;'' + 0(x5). 



(B.19) 
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